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ABSTRACT 


Maxwell's  equations  for  linear  media  are  reformulated  through 
linear  operator  and  generalized  transform  techniques  into  an  equivalent 
metric  integral  equation.  An  exp'" ;it  formal  solution  to  the  equation 
is  obtained  recursively,  providing  sequence  of  operations  to  be 
applied  to  the  electrical  parameters  of  the  medium  to  yield  the  charac¬ 
teristic  existence  conditions,  the  set  of  normal  modes,  and  the  electro¬ 
magnetic  fields  in  response  to  given  sources.  The  results  are  applicable 
to  time-invariant,  linear  media  which  may  be  inhomogeneous,  anisotropic, 
nonuniform,  dissipative,  dispersive,  with  any  source  distribution. 
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1.  INTRODUCTION 

Electromagnetic  phenomena  are  completely  described  mathemati¬ 
cally  by  Maxwell's  equations.  The  fundamental  problem  of  electromag¬ 
netic  theory  is  to  transform  the  implicit  description  embodied  in  the 
equations  into  explicit  formulations  of  these  phenomena  through  some 
procedure  for  solving  the  equations.  It  is  a  widely  recognized  truism 
that  Maxwell's  equations  are  amenable  to  solution  in  only  a  few  cases 
possessing  rather  special  simplifying  features.  A  less  restricted 
type  of  problem  is  usually  solved,  if  at  all,  by  appeal  to  methods 
often  specifically  designed  for  that  particular  problem.  There  remains 
a  wide  class  of  problems  that  resist  analysis.  A  general  method  of 
attacking  such  problems  is  developed  in  this  work. 

In  typical  problems  that  do  yield  to  standard  analysis,  the 
geometry  possesses  considerable  symmetry,  the  medium  consists  of  only 
a  few  homogeneous  regions,  the  coordinate  system  appropriate  to  the 
structure  of  the  system  is  one  of  the  few  in  which  the  wave  equation  is 
separable,  and  the  boundaries  are  of  such  nature  and  shape  that  it 
becomes  feasible  to  attempt  a  straightforward  solution  to  Maxwell's 
equations.  The  procedure  is  then  to  write  these  equations  and  the 
associated  boundary  conditions  for  each  region  of  the  space  and  to  solve 
each  such  set  of  partial  differential  equations  by  separation  of  vari¬ 
ables.  This  introduces  a  number  of  physically  meaningful  separation 
constants  and  arbitrary  coefficients.  The  individual  solutions  must 
then  be  made  to  fit  the  boundary  conditions,  particularly  those  arising 
at  the  interfaces  between  regions  of  differing  electrical  properties. 
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Hils  produces  an  auxiliary  set  of  equations,  no  longer  differential  but 
merely  algebraic  and  transcendental,  which,  beside  determining  the  arbi¬ 
trary  coefficients,  provides  a  relation  among  the  separation  constants. 
This  characteristic  relation,  which  may  be  a  dispersion  relation  for  a 
propagation  problem  or  a  resonance  condition  for  a  cavity,  is  often 
more  illuminating  and  useful  than  the  explicit  expressions  for  the  field 
components  themselves. 

Ihe  class  of  problems  for  which  solutions  in  closed  form  are 
obtainable  in  this  way  is  quite  small.  If  the  problem  is  not  sufficiently 
simple  to  permit  the  above  direct  attack,  recourse  is  often  had  to  some 
farm  of  perturbation  approach.  A  similar  but  simpler  problem  is  first 
solved  exactly  and  the  set  of  normal  modes  so  obtained  is  used  as  a 
basis  for  representing  the  solution  to  the  perturbed  case  as  a  superposi¬ 
tion  of  modes  of  the  simple  case.  Wie  extent  to  which  each  normal  mode 
contributes  to  the  expansion  remains  to  be  determined  through  the  use  of 
orthogonality  properties  enjoyed  by  these  modes  by  virtue  of  the  symmetry 
of  the  unperturbed  configuration.  Green's  functions,  the  use  of  which 
greatly  extends  the  class  of  problems  solvable  by  these  methods,  are 
usually  found  in  the  same  way. 

The  mathematical  elegance  associated  with  the  closed  forms  of 
solution  and  of  the  expressions  in  terms  of  more  or  less  well-known  and 
well-tabulated  functions  often  turns  out  to  be  an  illusory  advantage  of 
the  above  methods  of  solution.  The  characteristic  equations,  which 
provide  in  algebraic  form  implicit  information  about  physical  effects 
such  as  dispersion,  scattering,  or  resonance,  may  require  exhaustive 
studies  of  the  mathematical  properties  of  the  constituent  functions,  as 
is  also  the  case  when  the  distribution  and  flow  of  electromagnetic 
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energy  is  to  be  computed  from  the  field  components.  These  tasks,  though 
they  are  facilitated  by  a  host  of  recursion  and  orthogonality  properties 
possessed  by  these  functions,  must  often  be  relegated  to  automatic 
computers,  which  must  in  turn  rely  upon  the  power  series  expansions  of 
the  elementary  functions  involved. 

Abstracting  the  significant  features  of  these  general  and 
powerful  methods  of  solution,  the  following  observations  may  be  made. 
Firstly,  the  linearity  of  the  equations  underlies  the  methods,  heavy 
reliance  being  placed  upon  the  superposition  principle.  Secondly,  the 
search  for  closed  forms  of  solutions  may  be  self-defeating  in  that 
their  complexity  may  inhibit  their  interpretation  and  application. 
Thirdly,  the  recursive,  reciprocal,  and  orthogonality  properties  of  the 
solutions  so  obtained  are  actually  properties  of  the  normal  modes  asso¬ 
ciated  with  a  given  configuration  of  materials  in  the  medium,  not  just 
of  the  various  elementary  functions  involved,  which  are  after  all 
defined  as  solutions  of  the  separated  wave  equations.  Fourthly,  the 
algebraic  characteristic  equations  effectively  supplant  the  original 
differential  equations,  are  equivalent  to  the  Maxwell  equations  spe¬ 
cialized  to  the  particular  medium,  and  contain  essentially  the  same 
physical  information.  Finally,  these  characteristic  equations  are 
entirely  determined  by  the  electrical  constitution  of  the  medium  and 
should  be  obtainable  without  reference  to  the  associated  electromagnetic 
fields. 

These  principles  are  the  foundation  for  the  method  of  solution 


to  be  developed.  Primarily,  the  goal  is  to  reformulate  the  contents  of 
Maxwell's  equations  so  as  to  lead  from  a  description  of  the  spatial 
distribution  of  matter  in  the  medium  directly  to  a  single  characteristic 


k 


equation  which  prescribes  the  conditions  for  the  existence  of  a  mode* 
Secondarily,  a  formulation  is  sought  which  will  yield  the  conqplete  set 
of  normal  modes  for  any  linear  medium.  The  form  in  which  the  results 
are  to  be  expressed  is  to  be  left  arbitrary,  the  choice  to  be  made  in 
advance  on  the  basis  of  convenience  and  practicality. 

In  view  of  the  nondifferential  nature  of  the  characteristic 
equations  whose  content  is  equivalent  to  that  of  the  original  equations, 
the  present  approach  seeks  to  eliminate  the  differential  character  of 
the  equations  at  the  outset.  This  can  be  accomplished  by  means  of  a 
transformation  from  coordinate  space  to  some  new  domain  in  which  differ¬ 
entiation  is  replaced  by  other  operations  which  also  absorb  the  asso¬ 
ciated  boundary  conditions.  Pull  advantage  is  taken  of  the  linearity 
of  the  problem,  in  assuming  representations  of  unknown  functions  as 
general  summations,  in  expressing  the  response  to  several  excitations 
as  the  superposition  of  individual  responses,  and  as  the  essence  of  the 
reformulation  in  terms  of  linear  operator  techniques  from  which  general 
properties  of  electromagnetic  waves  may  be  deduced. 

A  guiding  principle  in  the  reformulation  of  Maxwell's  equations 
is  that  no  information  which  could  be  extracted  from  the  original  equa¬ 
tions  be  rendered  unobtainable  from  the  new  machinery.  This  requires 
that  all  information  contained  in  a  description  of  the  electrical  pro¬ 
perties  of  the  constituents  of  the  medium  be  retained  intact  through  the 
various  transformations,  operations,  and  manipulations  to  be  prescribed. 
In  effect,  there  will  be  presented  an  alternate  set  of  mathematical 
equations  to  that  of  Maxwell  through  which  the  electromagnetics  problem 
of  which  they  are  a  model  may  be  solved.  It  is  claimed  that  the  new 
set  of  equations  may  profitably  and  without  loss  be  considered  as  an 
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alternate  starting  point  for  the  investigation  of  a  large  class  of 
problems  in  electromagnetic  theory. 

Beyond  that,  it  will  be  shown  that  the  new  formulation  is 
subject  to  a  systematic  process  for  extracting  the  complete  solution  to 
the  problem.  Thus,  a  unified,  systematic,  formal  solution  to  Maxwell's 
equations  for  linear  media  will  be  presented.  The  solution  will,  of 
course,  be  only  formal,  for  the  machinery  developed  must  await  an  input 
in  the  form  of  a  precise  description  of  the  medium  for  which  the  solu¬ 
tion  is  desired  before  it  can  yield  the  dispersion  relation  or  other 
existence  conditions,  as  well  as  the  normal  modes  themselves,  as  the 
output. 

The  resolution  process  to  be  developed  may  become  exceedingly 
cumbersome  in  many  cases,  as  a  host  of  quadratures  may  be  called  for. 
This  will  be  a  manifestation  of  the  complex  nature  of  the  problems  to 
be  attacked  and  of  the  generality  of  the  method.  Automatic  computers 
of  a  high  degree  of  sophistication  may  well  be  required  to  render  the 
process  an  efficient  one.  However,  an  added  feature  of  the  method  is 
that  the  complexity  of  the  calculations  will  be  to  some  extent  under 
the  control  of  the  user. 

No  claim  is  made  to  adherence  to  strict  mathematical  rigor  in 
what  follows.  Hie  arguments  presented  are  intended  to  demonstrate  the 
plausibility  of  the  results.  Historical  precedents  for  various  methods 
to  be  employed  will  be  indicated  as  the  appropriate  stages  are  reached. 

A  computationally  trivial  illustration  of  the  theory  will  be  given  which 
is,  however,  sufficiently  general  to  be  of  interest.  The  example  will 
be  further  specialized  to  a  specific  medium  of  the  gyrotropic  type  in 
order  to  illustrate  the  nature  of  the  computations  required.  Finally, 
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possible  refinements  of  the  theory  will  he  indicated,  with  suggestions 
for  improving  and  extending  the  validity,  facility,  applicability,  and 
utility  of  this  reformulation  of  Maxwell's  equations. 
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2.  SCOPE  OF  PROBLEM 

The  problem  to  be  attacked  herein  is  that  of  determining  the 
characteristics  of  the  electromagnetic  fields  which  may  exist  in  a  given 
linear  medium  in  response  to  known  sources.  To  be  more  precise,  the 
restrictions  on  the  medium  are  that  there  be  a  definite,  linear,  time- 
invariant  relation  among  the  field  intensity,  flux  density,  and  current 
density  vectors  at  each  point  in  space.  These  conditions  are  sufficient¬ 
ly  weak  to  permit  consideration  of  a  medium  which  is  nonunifarm,  inhomo¬ 
geneous,  anisotropic,  lossy,  with  or  without  sources.  It  is  required 
merely  that  the  region  be  describable  electrically  by  giving  the  three 
tensor  fields  of  capacitivity,  permeability,  and  conductivity  at  every 
point.  The  problem  is  to  extract  frcm  this  data  the  characteristics  of 
the  macroscopic  electromagnetic  fields  which  may  be  excited  in  the  region 
of  interest  by  sources  within  or  without  it. 

The  time-invariance  imposed  upon  the  medium  permits  a  consid¬ 
erable  simplification  of  the  problem  in  that  a  harmonic,  steady-state 
analysis  will  suffice,  with  little  loss  of  generality.  In  addition, 
this  condition  permits  the  combination  of  two  of  the  three  tensor  fields, 
the  capacitivity  and  conductivity,  into  a  single  permittivity  tensor 
field.  Thus,  an  arbitrary  geometrical  configuration  of  lossy  dielectric 
and  permeable  materials  with  any  physical,  even  discontinuous,  variation 
of  electrical  properties  is  describable  by  two  dyadic  functions  of  posi¬ 
tion  and  frequency.  The  condition  that  the  relation  between  the  field 
intensities  and  flux  densities  be  an  unambiguous  one  imposes  one  final 
restriction  upon  the  systems  to  be  considered.  It  is  required  that  the 
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constitutive  dyadic  s,  both  the  permittivity  and  the  permeability,  have 
nonzero  eigenvalues. 

In  a  system  satisfying  the  conditions  outlined  above,  all  the 
relevant  electrical  information  concerning  the  system  will  be  contained 
in  the  two  constitutive  tensor  fields.  They  alone  are  sufficient  to 
prescribe  the  form  of  the  normal  modes  of  the  medium.  If  the  sources, 
both  the  electric  and  magnetic  current  densities,  be  specified  as  well, 
the  two  dyadic s  will  determine  the  actual  electromagnetic  field  config¬ 
uration  everywhere  in  space.  It  is  proposed  to  prescribe  a  set  of 
operations  for  transforming  the  information  contained  in  the  two  con¬ 
stitutive  dyadic  functions  into  information  about  the  electromagnetic 
fields  which  will  exist  in  the  given  medium. 

Under  the  stipulated  conditions,  the  macroscopic  electromag¬ 
netic  fields  and  source  currents  will  be  complex  vector  functions  of 
position,  time-harmonic  at  the  radian  frequency  ca,  and  describable  by 
three-component  vector  functions  in  some  arbitrarily  chosen  coordinate 
system.  Hie  electrical  constitution  of  the  medium  is  specified  by  three 
tensor  fields:  the  relative  capacitivity,  ^(r),  the  relative  permea¬ 
bility,  u(r),  and  the  conductivity,  o(r),  where  r  represents  the  posi¬ 
tion  vector.  By  virtue  of  the  nonuniformity  and  inhomogeneity  of  the 
medium,  these  tensors  are  functions  of  position  which  may  have  discon¬ 
tinuities,  for  example  at  the  interface  between  a  dielectric  and  free 
space.  If  the  medium  is  isotropic,  the  tensors  will  reduce  to  scalars. 

In  order  to  exhibit  Maxwell's  equations  explicitly,  the 
following  definitions  will  be  used.  They  are  appropriate  for  the 
rational  mks  system  of  units. 
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Permeability  of  free  space: 
Capacitivity  of  free  space:  tQ 

Speed  of  light  in  vacuo:  c 

Intrinsic  impedance  in  vacuo:  T) 

Wave  number  in  vacuo:  k 

Imaginary  unit:  J 


2  , 
c  “o*o  ■  1 

k  ■  <c/c 


Hie  permittivity  tensor  field  is  a  combination  of  the  relative  capaci¬ 
tivity  and  conductivity  tensors: 

e(r)  *  /f(r)  -  J ( 'H/k )cr ( r )  (2.1) 

Hie  permittivity  and  permeability  dyadic s,  t(r)  and  p(r),  have  been 
assumed  to  be  nonsingular. 

For  convenience,  the  electromagnetic  field  quantities  will  be 
expressed  in  the  following  form,  wherein  E(r),  M(r),  Ce(r),  and  Cm(r) 
are  complex,  vector  functions  of  position  and,  implicitly,  of  the  wave 
number. 


Electric  field 

E(r)ei<at 

(2.2) 

Magnetic  field 

(J/Wr)ej£at 

(2.3) 

Electric  current  density 

(k/jTl)«(r)Ce(r)eJ<Dt 

(2.h) 

Magnetic  current  density 

kp(r)Ca(r)eJtDt 

(2.5) 

In  terms  of  these  quantities 

,  Maxwell's  equations  take  the 

curl  E  »■  k  pM  - 

k“Cm 

(2.6) 

curl  M  -  k«E-k«Ce  (2.7) 

These  equations  are  to  be  satisfied  at  every  point  of  the  space.  Bound¬ 
ary  conditions,  other  than  the  regularity  of  the  solutions,  need  not  be 
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specified  since  they  have  been  incorporated  in  the  spatial  variation  of 
the  constitutive  tensors,  «(r)  and  n(r). 
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3*  OPERATOR  FORMULATION 

Maxwell's  equations  are  a  pair  of  coupled  vector  partial  dif¬ 
ferential  equations  with  variable  tensor  coefficients.  For  the  purpose 
of  exhibiting  a  general  solution  to  these  equations,  it  is  convenient 
to  reformulate  them  as  a  single  operator  equation.  This  eclipses  the 
multiple  and  coupled  character  of  the  equations  and  permits  considera¬ 
tion  of  Maxwell's  equations  as  an  abstract  relation  between  the  sources 
and  the  resultant  fields.  The  formulation  will  be  that  of  an  abstract 
operator  acting  upon  an  abstract  quantity  representing  the  response  to 
transform  it  into  one  that  represents  the  excitation.  Ultimately,  this 
will  facilitate  the  concretization  of  an  abstract  solution  to  this 
general  problem. 

Operator  methods  are  a  well-known  and  potent  tool  in  quantum 
mechanics,  but  less  common  in  electromagnetic  theory.  As  those  of  a 
linear  algebra,  such  methods  have  an  intrinsic  value  in  educing  certain 
invariant  properties  of  the  solutions  to  operator  equations,  such  as 
reciprocity^  or  the  reality  of  the  eigenvalues.  As  such,  these  operator 
methods  have  been  explored  by  Bresler,  Marcuvitz,  and  others  ”  for 
their  applicability  to  electromagnetic  theory  problems.  Herein,  these 
aspects  of  the  operator  calculus  will  not  be  pursued;  it  will  be  intro¬ 
duced  mainly  to  facilitate  the  formulation  of  an  abstract,  formal  solu¬ 
tion  to  Maxwell's  equations. 

From  the  operator  point  of  view,  Maxwell’s  equations  prescribe 
a  set  of  operations  to  be  performed  upon  the  electromagnetic  field  and 
source  vectors,  the  results  of  which,  when  equated,  express  the  physical 
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interrelation  among  the  fields.  As  the  fields  are  three-vectors,  the 
operators  acting  upon  them  should  be,  apart  from  scalar  multipliers, 
three-dyadies.  Hie  permittivity  and  permeability  tensors  are  such  dyad- 
ics;  when  applied  to  the  field  intensity  vectors,  they  produce  the  flux 
density  vectors.  Similarly,  insofar  as  it  produces  a  new  vector  field 
from  a  given  one,  the  curl  operation  should  be  expressible  as  a  dyadic 
operator. 


Hiis  is  indeed  possible,  in  any  particular  coordinate  system. 
Hius,  in  a  rectangular  coordinate  system,  (x,y,z),  the  three-dyadic 
representation  of  the  curl  operator  is 


curl 


a 

a 

0 

•  37 

SF 

a 

a 

0 

"  37 

a 

a 

3y 

a7 

0 

(3.1) 


In  a  cylindrical,  (p,< 6,z),  or  spherical,  (r,©,0),  coordinate  system, 

similar  dyadic  representations  are 

Cylindrical: 


curl 


0 

a 

-37 

i  a 

P 

a 

a 

37 

0 

"^p 

i  a 

i  a  „ 

7  ‘37 

p37p 

0 

(3-2) 
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Spherical: 


curl 


0 

1  a 

r  sin  6  9$ 

1  9 

”  r  "Se 


1  9 

r  sin  0 


0 

1  9 
r  "5? 


r 


1 

r  sin  0 
1 


^  sine 
9 


0 


(3.3) 


The  curl  operator  may  thus  he  thought  of  as  an  abstract  dyadic  operator 
which  may  be  given  an  explicit  expression  in  any  coordinate  system. 

There  may  now  be  introduced  two  immittance  operators  as  combi¬ 
nations  of  the  above  operators,  obtained  formally  by  the  rules  of  matrix 
inversion  and  multiplication.  The  impedance  operator  is  a  three-dyadic 
defined  as 

imp  «  (ks)_1curl  (3-M 

and  the  admittance  operator  is 

adm  »  (k(i)”*curl  (3*5) 

The  two  field  intensity  vectors  may  be  combined  into  a  single, 
abstract,  six-ccmponent  field  vector,  S(r),  and  the  two  excitation  vec¬ 
tors  into  a  source  six-vector,  C(r),  whose  structures  are  exhibited  in 
the  following  partitioned  column  vectors. 


rE(r)-| 

rcjrn 

s(r>  -L,J 

C(r)  - 

e 

The  Maxwell  operator  is  an  abstract  six-dyadic  operator  which 
will  operate  upon  such  six-vectors.  Its  structure  is  given  by  the 
following  partitioned  matrix. 


Ik 


maxi 


.  imp 


(3-7) 


Here,  1  represents  the  unit  three-dyadic. 

In  terms  of  the  above  abstract  quantities.  Maxwell's  equations 
become  the  single  operator  equation 

maxi  S(r)  -  C(r)  (3*8) 


The  six-dyadic,  differential  operator,  maxi,  incorporates  the  physical 
Interaction  among  the  electromagnetic  fields  as  well  as  the  electrical 
properties  of  the  medium. 

In  abstract  terms,  the  solution  to  the  problem  requires  the 
determination  of  the  operator  inverse  to  maxi,  that  is,  that  operator 
or  process  which,  when  applied  to  the  source,  C(r),  will  yield  the  fields, 
S(r). 
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k.  TRANSFORM  FORMULATION 

Hie  operator  form  of  Maxwell's  equations,  Eq.  (3.8),  involves 
a  differential  operator  acting  upon  the  unknown  function.  A  direct 
integration  of  this  differential  equation,  say  by  substitution  of  trial 
functions,  is  not  feasible  since  the  solution  cannot  be  expected  to  be 
expressible  in  closed  form  in  terms  of  known,  tabulated,  elementary 
functions.  Realization  of  this  raises  the  question  of  just  what  will 
constitute  a  solution  to  the  problem;  that  is,  how  may  the  six-compo¬ 
nent  vector  function  of  position,  S(r),  which  solves  the  equation  be 
exhibited? 

This  question  of  the  presentation  of  the  solution  may  be  re¬ 
solved  by  choosing  some  convenient  basis  of  representation  for  the 
function  to  be  described,  much  as  a  coordinate  system  is  chosen  for 
the  mensuration  of  the  space.  Hie  basis  will  consist  of  a  set  of  func¬ 
tions  with  the  property  of  completeness,  at  least  with  respect  to  the 
class  of  functions  which  are  sufficiently  regular  to  be  possible  solu¬ 
tions  to  the  physical  problem  here  considered.  The  unknown  solution 
will  then  be  expressed  as  an  expansion  in  this  complete  set.  After  the 
arbitrary  complete  set  has  been  selected,  the  problem  remaining  is  that 
of  determining  the  generalized  Fourier  coefficients  of  the  expansion. 

By  the  above  procedure,  the  problem  will  have  been  transformed 
from  the  solving  of  a  differential  operator  equation  in  position  space 
to  the  solving  of  a  corresponding  equation  in  the  space  of  the  Fourier 
coefficients.  The  spatial  variables  will  have  been  eliminated,  just  as 
the  Laplace  transform  typically  eliminates  the  time  variable  from  an 
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equation  by  transforming  to  the  frequency  domain.  Most  important,  the 
differential  operator  in  physical  space  will  have  been  replaced  by  one 
of  an  algebraic  or  integral  nature  in  transform  space. 

In  electromagnetic  theory,  a  form  of  the  process  contemplated 

Q 

here  was  developed  by  Schelkunoff'  to  produce  the  generalized  telegra¬ 
phist's  equations.  In  a  typical  problem  of  propagation  in  a  nonuniform 
waveguide,  expansion  functions  are  chosen  from  the  solution  of  a  related, 
simpler  problem,  say  that  of  a  limiting  case  of  the  actual  one.  The 
resulting  transformation  eliminates  the  transverse  variables,  so  that  a 
set  of  coupled  ordinary  differential  equations  in  only  the  axial,  variable 
is  obtained  for  the  coefficients  of  the  expansion.  These  telegraphist's 
equations  are  of  the  nature  of  transmission  line  equations  and  the  cou¬ 
pling  coefficients  provide  much  information  about  the  propagation  of 
modes  in  the  nonuniform  system.  This  process  is  somewhat  generalized 
herein,  where  the  expansion  functions  remain  arbitrary,  all  the  spatial 
variables  are  eliminated,  the  differential  character  of  the  equation  is 
entirely  suppressed,  and  vector  and  dyadic  functions  are  treated.  The 
procedure  will  yield  an  equation  for  the  transform  variables  to  which  a 
definite  set  of  operations  may  be  applied  to  extract  the  complete 
solution. 

®ie  translation  from  the  original  to  the  transform  space  must 
be  invertible;  that  is,  a  solution  obtained  in  transform  space  should  be 
subject  to  being  carried  back  to  the  original  domain,  if  desired.  There 
will  then  be  no  loss  of  information  attendant  to  the  consideration  of 
the  problem  in  the  transform  space. 

The  original  space  is  indexed  by  the  position  vector  r,  which 
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may  be  considered  a  composite  index  comprising  the  three  coordinates  of 
the  space.  In  the  transformation,  each  coordinate  or  spatial  index  is 
eliminated  at  the  cost  of  introducing  a  corresponding  transform  index, 
just  as  the  frequency  variable  supplants  the  time  variable  when  a  Fourier 
transform  is  applied  to  a  problem  in  the  time  domain.  The  transform 
space  will  thus  be  indexed  by  a  composite  index,  s.  The  transformation 
of  a  function  of  position,  g(r),  is  achieved  by  multiplying  it  by  the 
transformer  t(s,r)  and  integrating  throughout  the  space.  This  elimi¬ 
nates  the  position  index  r  and  leaves  the  transform  G(s),  a  function 
of  the  composite  index  s.  The  inverse  transformation  is  achieved  by 
multiplying  the  transform  G(s)  by  the  transformer  T(r,s)  and  inte¬ 
grating  or  summing,  as  appropriate,  over  each  variable  of  the  composite 
index,  s,  leaving  the  function  of  position,  g(r). 

Such  multiple  integrations  and  summations  as  are  involved  in 
the  above  processes  will  here  be  denoted  by  the  generic  summation  symbol 
$.  To  be  precise,  this  symbol  will  dictate  that  the  expression  which 
follows  it  be  integrated  or  summed  over  the  complete  range  of  the  vari¬ 
ables  which  are  contained  in  the  repeated,  dummy,  composite  indices 
appearing  in  the  summand.  Summation  will  be  implied  for  the  discrete 
variables  in  the  composite  index  and  integration,  possibly  with  appro¬ 
priate  weight  functions,  for  the  continuous  variables.  Urns,  the  above 
transform  pair  will  be  written 

G(s)  ■:  $  t(s,r)g(r)  g(r)  -  $  T(r,s)G(s)  (4.1) 

The  unit  function  for  such  summations  will  be  denoted  by  l(u,s).  It  has 
the  property  expressed  by 


G(u)  -  $  l(u,s)G(s) 


(4.2) 
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mis  idemfactor  thus  comprises  a  Dirac  delta  function  for  each  variable 
In  the  composite  Index  whose  range  is  a  continuum,  a  Kronecker  delta 
for  each  discrete  variable,  and  a  unit  dyadic,  as  required. 

As  an  Illustration  of  the  use  of  this  compact  notation,  con¬ 
sider  a  function  of  position,  g(r),  in  a  cylindrical  coordinate  system 
to  which  it  is  desired  to  apply  a  Fourier-Bessel  transform.  In  this 
case,  the  composite  index  r  represents  (p,<S, z)  and  the  transformation 
will  lead  to  a  space  indexed  by  a  »  (q,n,0).  lhe  transformer  t(s,r) 
and  Its  inverse  T(r,s)  are  given  by 

t(s,r)  -  t(q,n,0;p,<a,z)  *  e-d(n0+Pz)  Jn(qp)  (4.3) 

T(r,s)  -  T(p,$,z;q,n,0)  ■  (2*)“2  JQ(qp)  (4.4) 


The  ranges  of  p  and  q  are  0  to  ®,  with  weight  functions  p  and  q,  respec¬ 
tively;  the  range  of  0  is  0  to  2 *,  that  of  n  is  all  positive  and  negative 
integers;  the  ranges  of  z  and  0  are  -®  to  ®.  Accordingly,  Bqs.  (4.1) 

would  in  this  case  be  interpreted  as 
»  2«  ® 


O(q,n,0) 


g(P,0,*O 


Jdz  |d0  |pdp  t(q,n, 0; p,0,z)g(p,tf,z) 


|qdq  T(p,0,z;q,n,0)G(q,n,0) 


(4.5) 


For  the  problem  at  hand,  the  spatial  functions  and  their  trans¬ 
forms  are  six-vectors.  Correspondingly,  the  transformers  will  be  slx- 
dyadics.  In  fact,  let  there  be  chosen  some  convenient  dyadic  complete 
orthonormal  set  c(r,s)  with  its  inverse  d(s,r).  With  no  loss  of 
generality,  these  six-dyadlcs  may  be  taken  to  be  diagonal.  The  complete¬ 
ness  property  is  expressed  by 
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$  c(r,s)d(s,p)  >  l(r,p)  (4.6) 

and  that  of  orthonormality  by 

$  d(u,r)c(r,s)  -  l(u,s)  (4.7) 

The  transformation  law  for  a  column  vector  v(r)  and  the  inver¬ 
sion  of  its  transform  V(s)  will  then  be 

v(r)  o  $  c(r,s)V(s)  V(s)  =  $  d(s,r)v(r)  (4.8) 

The  transformation  of  a  row  vector  w(r)  is  to  be  performed  in  the  con¬ 
verse  manner. 

w(r)  -  $  W(s)d(s,r)  W(s)  -  $  w(r)c(r,s)  (4.9) 


Hie  properties  expressed  by  Eqs.  (4.6)  and  (4.7)  then  guarantee  that  the 
form  of  the  scalar  inner  product  of  two  vectors  will  be  preserved  by  the 
transformation . 

$  w(r)v(r)  a  $  W(s)v(s)  (4.10) 

Another  consequence  of  these  transformation  laws  is  that  a  dyadic  kernel 
operator  equation  in  position  space  such  as  x(r)  ■  $  A(r,p)v(p)  will  be 
translated  into  an  equation  of  similar  form  in  transform  space,  X(u)  ■ 
$  B(u,s)V(s),  provided  that  the  transformation  law  for  the  dyadic  oper¬ 
ator  A(r,p)  is  taken  as 

B(u,s)  »  $  d(u,r)A(r,p)c(p,s)  (4.11) 

It  follows  that  any  linear  dyadic  operator  in  position  space 
whose  nature  is  such  that  there  could  be  constructed  for  it  a  represen¬ 
tation  as  a  kernel  of  an  Integral  operator  will  have  a  representation  in 
transform  space  in  the  farm  of  a  dyadic  function  of  two  composite  trans¬ 
form  indices.  This  transform  of  the  operator  may  be  obtained  by  allowing 
the  position-space  operator  to  act  upon  c(r,s),  premultiplying  the 
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result  "by  the  dyadic  d(u,r)  and  integrating  throughout  the  space,  as 
in  Eg.  (4.11).  The  actual  kernel  representation  of  the  operator  in 
coordinate  space  need  not  he  found  explicitly;  indeed,  that  kernel  might 
he  highly  singular,  possibly  involving  delta  functions  and  their  deriva¬ 
tives.  This  situation  is  analogous  to  that  which  obtains  in  the  quantum 
mechanical  representation  of  dynamical  variables  as  operators 

The  hypothesis  is  now  made  that  the  abstract  maxi  operator  of 
Eq.  (3.8)  is  such  an  operator.  The  rigorous  justification  of  this  hy¬ 
pothesis  will  be  held  in  abeyance  while  its  consequences  are  explored. 

Proceeding  now  with  the  transformation  of  the  operator  form  of 
Maxwell's  equations,  let  P(s)  and  Q(s)  be  the  transforms  of  the  vectors 
S(r)  and  C(r),  respectively.  That  is. 


S(r) 

-  $  c(r,s)F(s) 

F(s) 

■  $  d(s,r)S(r) 

(4.12) 

C(r) 

-  $  c(r,s)Q(s) 

Q(s) 

=  $  d(s,r)C(r) 

(4.13) 

In  accordance  with  the  preceding  discussion,  let  the  transform  of  the 
maxi  operator  be  D(u,s). 

D(u,s)  -  $  d(u,r)  maxi  c(r,s)  (4.14) 

Upon  premultiplying  both  sides  of  Eq.  (3*8)  by  d(u,r),  substituting  the 
expansion  for  S(r),  and  integrating.  Maxwell's  equations  become  the 
Operator  transform  equation 

$  D(u,s)F(s)  -  Q(u) 


(4.15) 
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5.  SIGNIFICANCE  OF  REFORMULATION 

The  equation  obtained  here  for  the  transforms  of  the  electro¬ 
magnetic  fields,  Bq.  (4.15), 

$  D(u,s)F(s)  ■  Q(u) 

merits  closer  examination.  Despite  the  notatlonal  disguise,  it  is 
entirely  equivalent  to  the  original  Maxwell  equations. 

In  this  equation,  F(s)  is  an  unknown  six-vector,  the  transform 
of  the  desired  field  vector  S(r).  The  six-vector  Q(s)  is  the  known 
transform  of  the  given  source  vector  C(r).  D(u,s)  is  a  six-dyadic 
kernel,  a  function  of  two  composite  indices  of  the  transform  space  and, 
implicitly,  of  the  wave  number,  k,  as  a  parameter.  It  Incorporates  in 
an  intimate  combination 

1)  the  physical  law  of  interaction  of  electromagnetic  fields,  namely 
that,  in  the  steady  state,  one  field  determines  the  spatial  rate 
of  change  of  the  other,  as  expressed  by  the  curl  operation; 

2)  the  electrical  constitution  of  the  medium  filling  the  space,  as 
described  by  the  tensor  capacitivity,  permeability,  and  conductiv¬ 
ity  at  every  point; 

3)  the  geometrical  configuration  of  the  material  bodies  whose  arrange¬ 
ment  in  space  is  incorporated  in  the  spatial  variation,  particularly 
the  discontinuities,  of  the  constitutive  tensors; 

k)  the  boundary  conditions  which  would  have  to  be  appended  to  the  orig¬ 
inal  partial  differential  equations  to  prescribe  the  field  disconti¬ 
nuities  at  Interfaces  between  regions  of  different  electrical 


properties; 
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5)  the  basis  of  representation  selected  to  egress  the  spatial  varia¬ 
tion  of  the  fields,  as  embodied  in  the  dyadic  transformers  c(r,s) 
and  d(s,r). 

Ike  summation  operation  gives  the  equation  the  character  of  a  generalized 
integral  equation  of  the  first  kind,  with  Q(s)  as  the  forcing  function, 
F(s)  as  the  unknown  function,  and  D(u,s)  as  the  kernel. 

A  review  of  the  procedure  involved  in  obtaining  the  dyadic  ker¬ 
nel  D(u,s)  will  show  that  no  information  about  the  electromagnetic  system 
has  been  lost  in  the  transition  from  the  original  form  of  Maxwell's  equa¬ 
tions  to  the  present  formulation.  An  appropriate  coordinate  system  for 
the  space  having  been  chosm,  the  medium  is  first  described  mathematically 
by  giving  the  two  constitutive  dyadics,  e(r)  and  n(r),  at  each  point  in 
space.  These  two  constitutive  matrices  are  inverted  and  then  premulti¬ 
plied  into  the  dyadic  representative  of  the  curl  operator  to  form  the 
immittance  operators  imp  and  adm,  as  in  Eqs.  (3.4)  and  (3.5).  These 
three-dyadic  operators  are  arranged  in  a  six-dyadic  as  in  Eq.  (3.7)  to 
create  the  maxi  operator  for  the  medium  in  question.  A  diagonal  dyadic 
complete  set  c(r,s)  and  its  inverse  d(s,r)  are  selected  as  a  basis  of 
representation.  The  maxi  operator  is  applied  to  c(r,s)  and,  finally, 
the  integrations  prescribed  in  Eq.  (4.14)  are  performed  to  produce  the 
kernel  D(u,s).  As  each  of  these  steps,  in  particular  the  transformation, 
is  reversible,  it  is  clear  that  all  the  pertinent  information  about  the 
system  has  been  retained  in  the  process  of  obtaining  D(u,s).  The  problem 
may  hence  be  solved  Just  as  thoroughly  by  an  attack  on  the  new  equation, 
Eq.  (4.13),  as  on  the  original  Maxwell  equations. 

The  most  important  characteristic  of  D(u,  s)  as  a  replacement 
far  the  maxi  operator  is  the  absence  of  spatial  coordinates  and  of 
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differential  operators.  A  major  advantage  of  the  reformulation  of 
Maxwell's  equations  lies  in  that  as  expressed  in  this  integral  equation 
form,  approximation  methods,  variational  techniques,  perturbation 
approaches,  and  abstract  operator  analyses  may  be  applied  to  the  problem 
more  easily  than  in  the  original  formulation  as  coupled,  vector,  partial 
differential,  variable  coefficient  equations.  In  fact,  a  complete,  exact 
solution  to  the  transformed  equation  may  be  extracted  through  the  appli¬ 
cation  of  a  definite  set  of  operations  to  the  known  kernel  D(u,s)  and 
source  Q(s).  Any  solution,  however  obtained,  for  F(s)  constitutes  the 
expansion  coefficients  for  the  field  vector  S(r)  in  the  selected  basis 
of  representation  c(r,s).  The  solution  F(s)  may  be  transformed  back  to 
the  desired  S(r),  in  accordance  with  Eq.  (4.12),  to  obtain  both  the 
electric  and  magnetic  fields  at  every  point  in  space. 

This  last  step  of  transformation  back  to  physical  space  will 
often  be  superfluous  as  much  information  may  be  gleaned  from  the  trans¬ 
form  itself.  For  example,  the  existence  of  transverse  electric  or  trans¬ 
verse  magnetic  waves  may  be  deduced  from  a  simple  examination  of  the 
transform  F(s)  for  any  vanishing  components.  Multiple  solutions  at  a 
given  frequency  are  likewise  reflected  in  the  multiplicity  of  the  solu¬ 
tion  for  the  transform  F(s)  at  that  wave  number. 

Less  trivially,  in  the  case  of  a  sourceless  region  for  which 
the  normal  modes  are  desired.  Maxwell's  equations  in  both  the  original 
and  transformed  forms  will  be  homogeneous  and  solutions  will  exist  only 
under  certain  conditions.  These  conditions  will  apply  Just  as  well  for 
the  solution  of  the  homogeneous  equation  in  transform  space,  typically, 
these  existence  conditions  determine  the  eigenfrequencies  of  the  system 
by  liqposlng  a  restriction  on  the  wave  number,  k.  Bjus  if  the  ay stem 
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tinder  consideration  is  a  cavity,  its  resonant  frequencies  will  be  ob¬ 
tained  in  solving  the  transformed  equations,  without  the  necessity  of 
transforming  back  to  coordinate  space.  Similarly,  if  propagation  of 
normal  modes  on  a  guiding  structure  is  considered,  the  conditions  for 
the  existence  of  solutions  to  the  homogeneous  transformed  equation  will 
constitute  the  dispersion  relation  for  the  medium.  In  a  propagation 
problem  the  propagation  constant  appears  simply  as  another  parameter  and 
the  characteristic  existence  condition  will  prescribe  a  relation  between 
the  wave  number  and  the  propagation  constant. 

Still  more  Information  about  the  system  may  be  extracted  from 
the  transformed  solution  by  a  limited  form  of  Inversion  of  the  trans¬ 
forms.  Ihus,  certain  quadratic  forms  in  the  fields  may  be  evaluated  in 
terms  of  the  transforms  without  first  obtaining  the  fields  explicitly. 
Same  quadratic  forms  of  importance  are 


We(r)  -  E+(r)«(r)E(r) 

(5.1) 

Wm(r)  «  M+(r)u(r)M(r) 

(5.2) 

N(r)  ■  E*(r)  x  M(r) 

(5.3) 

where  *  Indicates  the  complex  conjugate  and  +  the  Hermltlan  conjugate  of 
a  quantity.  The  Interest  of  these  quantities  arises  from  the  relations 
time-average  complex  electric  energy  density  -  («0/U)Wfi 
time-average  complex  magnetic  energy  density  ■  (co/4)Wm 
time-average  complex  Poynting  vector  ■  (l/2jT|)N* 

In  order  to  calculate  these  expressions  with  only  a  partial  inversion  of 
the  transforms  back  to  position  space,  some  of  the  dyadic s  to  be  used 
will  be  partitioned,  in  the  manner  of  Eqs.  (5.6), as  follows. 
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the  three  quadratic  forma  of  interest  will  be  obtainable  directly  from 
the  transforms  of  the  fields  througi 

V@(r)  -  *V+(u)We(u,r,s)V(s)  (5.12) 

Wn(r)  .  *  I+(u)w“(u,r,s)l(s)  (5.13) 

H(r)*n(r)  -  $  V+(u)Pn(u,r,s)l(s)  (5.14) 

It  may  be  noted  from  Eqs.  (5.12)  and  (5.I3)  that  a  passive 
medium  vill  be  characterized  by  dyadic s  W6  and  W®  that  are  positive  defi¬ 
nite  and  that  this  property  is  directly  related  to  the  positive  definite¬ 
ness  of  the  constitutive  dyadic s,  as  may  be  seen  from  Bqs.  (5.9)  and 
(5.10). 

Less  trivial  qualitative  and  semi -quantitative  information 
about  the  system  vill  be  obtainable  from  an  examination  of  the  character 
of  the  Poynting  vector  field.  This  may  be  inferred  from  Bq.  (5.14),  in 
which  n  will  typically  be  a  unit  vector  in  some  direction  of  Interest. 
Thus,  a  Poynting  vector  field  which  is  found  to  have  only  an  axial  compo¬ 
nent  will  evidently  characterize  a  propagating  wave.  Similar  qualitative 
characterizations  may  be  made  of  such  phenomena  as  reflection,  scattering, 
radiation,  backward  waves. 

Thus,  the  transform  equation,  Bq.  (4.15)>  is,  in  principle, 
entirely  equivalent  in  content  to  the  original  Maxwell  equations.  It 
has  the  merit  that  its  solution  may  be  systematized. 


6.  INTEGRAL  EQUATION  FORMULATION 


As  pointed  out  earlier,  the  operator  transform  equation  for 
the  electromagnetic  fields  obtained  herein  has  the  aspect  of  a  somewhat 
generalized  integral  equation  of  the  first  kind.  Such  equations  are 
generally  intractable  by  their  very  nature .  Furthermore,  the  quantities 
involved  in  the  equation  are  dyadics  of  sixth  order.  In  effect,  then, 
the  transform  equation  in  its  present  form,  Eq.  (4.15),  may  fairly  be 
said  to  be  unmanageable. 

The  six-dyadic  kernel  D(u,s)  has,  however,  been  so  constructed 
that  a  considerable  simplification  of  the  equation  is  possible.  The 
order  of  the  dyadics  to  be  manipulated  may  firstly  be  reduced  from  six 
to  three,  at  the  cost  of  introducing  a  pair  of  coupled  equations  for  the 
individual  field  transforms.  But  then  these  equations  may  be  recombined 
into  one  equation,  still  with  third-order  dyadics,  of  the  form  of  an  inte 
gral  equation  of  the  second  kind,  which  is  amenable  to  analysis.  There 
will  then  remain  merely  to  adapt  the  classical  Fredholm  theory  for  such 
equations  to  the  somewhat  generalized  type  of  Integral  equation  to  be 
obtained  here.  This  process  will  lead  ultimately  to  a  complete,  formal 
solution  to  Maxwell's  equations. 

The  reformulation  of  the  equations  continues  now  with  the 
further  dissection  of  the  six-dyadics  into  partitioned  matrices.  To 
accompany  the  partition  of  the  field  fector  F(s)  in  Bq.  (5«1»),  the 
source  vector  Q(s)  partitions  as 
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RKsH 

Q(s)  -  (6.1 

LJ(s)J 

In  view  of  the  construction  of  the  maxi  operator,  Eq.  (3*7),  and  of 
the  orthonormality  of  the  transformers,  e(r,s)  and  d(s,r),  as  expressed 
by  Eq.  (4.7),  the  six-dyadic  kernel  D(u,s)  has  the  following  highly 
significant  structure. 

D(u,s)  ■ 

The  three-dyadic  lnraittance  kernels,  Z(u,s)  and  Y(u,s),  are  defined  hy 
this  equation. 

Introducing  these  partitioned  farms  of  the  matrices  into  the 
transform  equation,  Eq.  (4.15),  it  is  seen  that  the  equation  can  he 
decomposed  into  the  network  equations 

V(u)  -  W(u)  +  $  Z(u,s)l(s)  (6.3) 

I(u)  «  J(u)  +  $  Y(u,s)V(s)  (6.4) 

In  these  coupled  equations,  the  transformed  field  three-vectors  V(s)  and 
I(s)  are  the  unknowns;  the  three-vectors  W(s)  and  J(s)  are  transforms  of 
given  sources;  the  Impedance  and  admittance  kernels  Z(u,s)  and  Y(u,s)  are 
known  three-dyadics,  obtained  by  inspection  of  Eq.  (6.2),  A  comparison 

Is  Invited  of  these  equations  with  the  static  network  equations  of 

14  Q 

Marcuvitz  or  with  the  dynamic  telegraphist's  equations  of  Bchelkunoff  . 

The  final  phase  of  the  reformulation  process  will  now  be 

expounded,  with  a  view  to  uncoupling  the  two  equations  and  producing  an 

equation  much  like  an  Integral  equation  of  the  second  kind,  for  which  a 


l(u,s) 

-  Y(u,s) 


-  Z(u,s) 
l(u,s) 


(6.2) 
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general  theory  may  he  formulated.  The  uncoupling  process  is  analogous 
to  that  which,  for  a  homogeneous  medium,  produces  the  wave  equation  from 
the  coupled  Maxvell  equations;  it  involves  simply  the  substitution  of 
Bq.  (6.4)  into  Eq.  (6.3). 

Some  preliminary  definitions  will  be  found  to  be  convenient. 

It  is  profitable  to  Introduce  a  new  parameter  into  the  equations,  for 
ease  of  analysis.  While  this  may  be  done  entirely  arbitrarily,  it  is 
particularly  convenient,  in  case  the  constitutive  tensors  are  independ¬ 
ent  of  frequency,  to  let  this  new  parameter,  X,  absorb  and  fill  the  role 
of  the  wave  number,  k.  This  will  facilitate  the  physical  interpretation 
of  some  of  the  equations  to  be  obtained.  It  should  be  emphasized  that 
the  present  restriction  to  frequency-independent  media  is  made  merely 
far  convenience  of  interpretation  and  will  be  lifted  later. 

The  system  being  studied  possesses  some  characteristic  dimen¬ 
sions.  Prom  these,  let  there  be  selected  a  convenient  quantity.  A, 
with  the  dimensions  of  an  area,  to  characterize  the  physical  size  of  the 
system.  The  dimensionless  parameter,  X,  is  defined  as 

\  m  l/k2A  (6.5) 

The  three-dyadic  kernel,  K(u,s),  is  defined  by 

K(u,s)  -  k2A  $  Z(u,v)Y(v,s)  (6.6) 

A  review  of  the  steps  Involved  In  calculating  the  two  inndttanee  kernels 
Z(u,s)  and  Y(u,s)  will  readily  show  that,  in  the  case  of  a  frequency- 
independent  medium,  the  kernel  K(u,s)  Is  independent  of  the  wave  number 
k,  and  hence  of  the  parameter  X.  Finally,  the  compound  source,  U(s),  is 
given  by 


U(u)  -  W(u)  +  *  Z(u,s)J(e) 


(6.7) 


A  direct  substitution  of  Eg.  (6.4)  into  Ec[.  (6.3)  new  yields 
the  operator  transform  reformulation  of  the  electromagnetics  problem  in 
generalised  Integral  equation  form, 

V(u)  -  U(u)  +  \  $  K(u,s)V(s)  (6.8) 

This  equation  is  to  be  solved  for  the  unknown  three-vector  V(s)  for  a 
given  compound  source  U(s),  kernel  K(u,s)  and  parameter  X.  When  the 
solution  is  introduced  into  Bj.  (6.4)  to  calculate  1(a),  the  complete 
solution  for  the  transformed  fields  is  obtained. 
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7.  METHOD  OP  SOLUTION 

(Hie  final  reformulation  of  Maxwell's  equations  obtained  here, 

»1.  (6.8), 

V(u)  -  U(u)  +  X.  $  K(u,s)V(s) 

has  the  aspect  of  a  linear  Fredholm  Integral  equation  of  the  second  kind. 
The  classical  Fredholm  theory  for  such  equations  must,  however,  be  modi¬ 
fied  for  the  case  at  hand  since  the  equation  is  somewhat  more  general 
than  the  prototype  equation  treated  by  Fredholm.  In  the  present  case, 
a  generalized,  multiple  summation  replaces  the  single  integration  of  the 
original  Fredholm  equation,  the  range  of  summation  is  generally  not  a 
simple  closed  Interval  of  integration,  and  dyadic  quantities  appear  in 
the  summands. 

This  last  innovation,  represented  by  the  dyadic  character  of 

the  kernel,  requires  a  nontrivial  modification  of  the  Fredholm  theory. 

For  tills  theory  ?  is  usually  developed  in  terms  of  determinants  with 

elements  formed  from  the  values  of  the  kernel  at  various  points.  Were 

this  process  carried  out  formally  in  the  present  case,  there  would  be 

obtained  determinants  with  nonccmmiting,  metric  elements,  whereupon  the 

theory  would  break  down.  A  formulation  is  required  which  avoids  the 

formation  of  determinants  and  is  applicable  to  this  metric  equation. 

Such  a  formulation  can  be  developed,  in  close  analogy  to  that 
17 

presented  by  Smithies,  1  but  with  certain  modifications  demanded  by  the 
dyadic  character  of  the  kernel.  The  quantities  normally  given  by  deter¬ 
minants  will,  in  this  metric  case,  be  expressed  in  terms  of  a  set  at 
recursion  relations. 
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I&e  method,  of  solution  to  be  presented  here  takes  full  advan¬ 
tage  of  the  linearity  of  the  equation  by  providing  a  solution  for  the 
case  of  just  a  unit  source  function.  The  solution  for  an  arbitrary 
source  can  then  be  obtained  by  superposition.  The  resolution  of  the 
problem  posed  by  Eq.  (6.8)  reduces  then  to  the  search  for  the  resolvent 
kernel  H(u,s;X),  which  plays  the  role  of  a  Green's  function  for  the 
equation,  in  terms  of  which  the  solution  to  the  equation 

V(u)  -  U(u)  +  X  $  K(u,s)V(s)  (7.1) 

far  an  arbitrary  source  U(s)  will  be 

V(u)  -  U(u)  +  X  $  H(u,s;X)U(s)  (7.2) 

Upon  substituting  this  assumed  solution  into  Eq.  (7.1),  noting  that  it  is 
to  be  satisfied  for  all  sources  U(s),  it  is  found  that  H(u,s;X)  must 
satisfy  the  resolvent  equation 

H(u,s;X)  *  K(u,s)  +  X  $  K(u,v)H(v,s;X)  (7.3) 

The  resolvent  H(u,s;X)  is  a  three-dyadic  kernel  associated 
with,  and  entirely  determined  by,  the  kernel  K(u,s)  of  the  equation. 

A  study  of  its  properties  is  equivalent  to  an  analysis  of  the  original 
electromagnetics  problem.  In  fact,  solving  Eq.  (7.3)  for  H(u,s;X)  will 
yield  the  solutions  for  the  transformed  fields  successively  flrcm  Eqs. 
(7.2)  and  (6.4).  Some  of  the  more  significant  properties  of  the  resol¬ 
vent  are  expressed  by  the  following  relations,  which  are  fairly  readily 
derivable  from  Eq.  (7.3) . 

$  K(u,v)H(v,s;X)  m  $  H(u,v;X)K(v,s) 

H(u,sjO)  =  K(u,s) 

-  $  H(u,v;X)H(v,s;X) 
ok 


(7.4) 

(7.5) 

(7.6) 
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H(u,sjk)  -  H(u,sj§)  -  (X-5)  $  H(u,v;X)H(v,s;S)  (7*7) 

In  Bq.  (7.6),  no  summation  over  X  is  Intended,  of  course,  although  it  is 
repeated  in  the  summand;  it  is  merely  a  parameter,  not  a  dummy  composite 
index. 

Another  important  property  of  the  resolvent,  following  immedi¬ 
ately  from  Eq.  (7.2),  is  that  as  long  as  H(u,s;X)  is  finite,  the  homo¬ 
geneous  equation 

V(u)  -  X  $  K(u,s)V(s)  (7.8) 

will  have  only  the  trivial  solution  V(s)  «  0.  The  homogeneous  equa¬ 
tion  will  have  nontrivial  solutions,  however,  but  only  for  certain 
characteristic  values  of  X.  At  these  values,  H(u,s;X)  will  have  to 
become  infinite. 

In  order  to  gain  greater  insight  into  the  nature,  structure, 
and  properties  of  the  resolvent,  it  is  instructive  to  consider  a  very 
special  but  important  case  of  Bq.  (7.1).  Ibis  is  that  of  a  degenerate 
kernel,  in  which  case  it  is  possible  to  give  an  explicit  expression,  in 
closed  farm,  for  the  resolvent.  A  degenerate  kernel  is  one  that  can  be 
factored  so  as  to  separate  its  two  composite  indices.  Explicitly,  an 
n  x  n  dyadic  kernel  K(u,s)  is  degenerate  if  it  can  be  expressed  in 
factored  form  as 

K(u,s)  -  A(u)B(s)  (7.9) 

where  A(s)  is  an  n  x  r  dyadic  function  of  just  one  composite  index  and 
B(s)  is  an  r  x  n  dyadic  function  of  one  index.  Let  the  obverse  of  the 
degenerate  kernel  be  defined  as  the  r  x  r  constant  dyadic 

R  -  $  B(s)a(s) 


(7.10) 


In  the  degenerate  case,  then,  Eq.  (7.1)  becomes 

V(u)  -  U(u)  +  \  A(u)  $  B(s)V(s)  (7.H) 

from  which  there  may  be  obtained  a  simple  metric  equation  for  the  auxil¬ 
iary  unknown  X  ■  $  B(s)V(s)  by  premultiplying  by  B(u)  and  summing: 

X  »  $  B(u)U(u)  URX  (7.12) 

This  equation  may  immediately  be  solved  for  X,  which  may  in  turn  be  sub¬ 
stituted  in  Eq.  (7»H)  to  obtain  the  solution  for  V(s).  Upon  comparing 
the  result  with  Eq.  (7.2),  it  is  seen  that  the  resolvent  for  this  degen¬ 
erate  case  is  given  by 

H(u,s;X.)  =  A(u)(l  -  \  R)"1  B(s)  (7.I3) 

The  structure  of  this  explicit  expression  for  the  resolvent  for 
a  degenerate  kernel  clarifies  many  of  the  properties  of  resolvents  for 
the  general  case.  The  properties  stated  in  Eqs.  (7.k)  -  (7.7)  may  be 
readily  verified  for  this  resolvent.  Of  greater  importance  for  the 
sequel  is  the  fact  that  the  resolvent  contains  as  a  common  denominator 
the  determinant  det  (1  -  \  R).  This  scalar  function  of  the  parameter 
X  determines  the  poles  of  the  resolvent  considered  as  a  function  of  \. 

The  secular  equation  is 

det  (1  -  \  R)  »  0  (7.1U) 

and  its  roots  are  those  values  of  X.  for  which  the  homogeneous  equation, 
Eq.  (7*8),  has  nontrivial  solutions.  In  fact,  it  may  be  seen  from  Eqs. 
(7.12)  and  (7.U)  that  the  solutions  to  the  homogeneous  equation  will  be 

V(s)  -  A(s)  X  (7.I5) 

where  X  is  any  eigenvector  of  the  obverse  matrix,  R. 

The  secular  equation,  Eq.  (7.1k),  is  a  polynomial  equation  of 
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degree  r,  so  that  there  will  he  r  characteristic  values  of  X,  count¬ 
ing  multiplicities.  This  assumes  that  the  order,  r,  of  the  obverse 
matrix  is  also  its  rank;  if  R  is  singular,  there  will  he  fever  roots, 
hut  then  a  different  factorization  of  the  kernel  could  have  been  found 
far  which  the  obverse  matrix  would  have  had  a  lower  order. 

These  results  for  the  degenerate  case  will  he  generalized  for 
that  of  a  general,  nondegenerate  kernel.  The  central  result  will  be,  as 
may  already  be  anticipated,  that  the  characteristic  values  of  X  will  he 
given  by  the  roots,  not  of  a  polynomial,  but  of  an  infinite  power  series 
in  X.  The  explicit  expression  for  the  resolvent  in  the  degenerate  case 
given  in  Eq.  (7.13)  provides  the  clue  to  the  structure  of  the  general 


resolvent . 
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8.  STRUCTURE  OP  RESOLVENT 


Eqs  (7«13)  is  a  prescription  for  calculating  the  resolvent  of 
a  degenerate  kernel.  The  poles  of  the  resolvent,  which  correspond  to 
nontrivial  solutions  of  the  homogeneous  equation,  occur  at  the  roots  of 
the  determinant  det(l  -  X  R),  a  polynomial  of  degree  r  in  X.  These 
characteristic  values  of  X  are  an  intrinsic,  though  latent,  property  of 
the  kernel  and  the  fact  that  there  is  a  finite  number,  r,  of  such  roots, 
counting  multiplicities,  is  the  essential  characteristic  of  a  degenerate 
kernel. 

A  nondegenerate  kernel  will  inherently  possess  an  infinite 
number  of  characteristic  values  of  the  parameter  X.  No  factorization 
into  matrices  of  finite  order  as  in  Bn.  (7»9)  will  then  be  possible.  It 
is  always  possible,  however,  to  approximate  the  kernel  arbitrarily  close¬ 
ly  in  such  factored  form,  provided  matrices  of  infinite  order  are  admit¬ 
ted.  This  is  so  because  this  decomposition  will  then  be  nothing  mare 
than  an  infinite  sum  of  products  of  functions  of  each  composite  index 
separately,  which  certainly  can  represent  any  sufficiently  regular  ker¬ 
nel  with  any  degree  of  accuracy.  The  obverse  matrix  of  a  nondegenerate 
kernel  so  expressed  will  be  of  infinite  order,  the  determinant  det(l-XR) 
will  be  not  a  polynomial  but  an  infinite  power  series  in  X,  and  the  num¬ 
ber  of  characteristic  roots  will,  consistently,  be  infinite.  The  nonde¬ 
generate  case  can  thus  be  considered  a  limiting  case  of  that  of  a  degen¬ 
erate  kernel,  but  with  the  transition  of  the  order,  r,  to  infinity. 

Evidently,  however,  the  calculation  of  the  resolvent  of  a  non¬ 
degenerate  kernel  by  Eq.  (7.13)  with  the  obverse  matrix,  R,  and  the  two 


factors,  A(u)  and  B(s),  of  infinite  order  is  not  feasible.  Nevertheless, 
the  fact  that  the  resolvent  could  he  so  expressed,  in  a  formal  sense, 
provides  the  clue  to  the  solution  of  the  problem  for  general  kernels. 

Far  this  possibility  of  considering  a  nondegenerate  kernel  as  a  limiting 
case  of  a  degenerate  one  Indicates  that  the  structure  of  the  resolvent 
will  be  that  of  a  resolvent  for  a  degenerate  kernel  and,  in  fact,  any 
property  of  a  resolvent  for  a  degenerate  kernel  which  does  not  depend 
explicitly  on  the  order  r  or  on  its  finiteness  will  be  possessed  by 
the  resolvent  for  the  general  kernel.  Recognition  of  this  fact  leads  to 
the  solution  of  the  problem. 

Accordingly,  in  direct  analogy  with  the  degenerate  case,  the 
resolvent  for  a  general  kernel  will  be  assigned  the  structure 

H(u,sj\)  -  C(u,sj\)/p(\)  (8.1) 

The  determinator,  p(\),  is  a  scalar  function  of  only  the  parameter  \  and 
is  the  analog  of  the  determinant  det(l  -  \  R)  of  the  degenerate  case. 

The  character! zer,  C(u,s;\),  is  a  three-dyadic  function  of  two  composite 
indices  and  of  it  corresponds  to  the  matrix  A(u)Q(k)B(s)  of  the 
degenerate  case,  where  Q(\)  is  the  adjoint  of  the  matrix  1  -  \  R. 

Both  the  characterizer  and  the  determinator  are  unknown.  Their 
relation  to  the  kernel  is  readily  determined  frcm  the  resolvent  equation, 
Bq.  (7.3),  by  multiplying  by  p(\).  The  resulting  characterizer  equation 
is 

C(u,sjk)  -  K(u,s)p(\)  +  \  $  K(u,v)C(v,s;\)  (8.2) 

This  is,  of  course,  insufficient  to  determine  both  C(u,s;\)  and  p(\). 

A  further,  though  still  partial,  specification  of  the  determinator  may 
be  made  in  analogy  with  the  degenerate  case  by  fixing  the  arbitrary 
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multiplicative  constants  of  C(u,s;\)  and  p(x.)  in  accordance  vith 

P(0)  -  1  (8.3) 

It  follows  fro©  Eq.  (7.5)  that 

C(u,sjO)  »  K(u,s)  (8.k) 

Further  progress  toward  the  determination  of  C(u,sj\)  and  p(\) 
may  he  made  by  prescribing  the  structure  of  the  determinator  to  be  the 
analog  of  that  of  the  degenerate  case;  that  is,  p(X)  is  to  be  a  polynomial 
in  \,  but  of  infinite  degree: 

p(>0  ’  Jo  V-”  (8.5) 

Eq.  (8.3)  prescribes  that  pQ  -  1;  all  other  coefficients  remain  unknown. 
Similarly,  the  character! zer  will  be  expanded  as 

C(u,s;\)  *  Cn(u,s)\  (8.6) 

Here,  CQ(u,s)  »  K(u,s)  and  the  other  three-dyadic  coefficients  are  un¬ 
known.  Upon  substitution  of  these  assumed  expansions  into  the  charac- 
terlzer  equation,  Eq.  (8.2),  there  is  obtained  the  recursive  character¬ 
izer  equation 

Cn(u>s)  -  K(u,s)pn  +  $  K(u,v)Cn-1(v,s)  (8.7) 

Eq.  (8.7)  is  a  recursion  relation  for  the  characterizer  coeffi¬ 
cients,  but  it  requires  a  knowledge  of  the  determinator  coefficients,  pQ, 
for  the  successive  calculation  of  the  Cn(u,s),  starting  from  the  known 
initial  coefficients 

CQ(u,s)  -  K(u,s)  P0  -  1  (8.8) 

There  remains,  then,  to  specify  the  determinator  coefficients,  whereupon 
all  the  characterizer  coefficients  will  be  obtainable  in  succession  from 
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the  recursive  characterizer  equation,  and  from  these  the  complete  solu¬ 
tion  for  the  general  kernel.  But,  at  this  point,  an  impasse  appears  to 
have  been  reached  since  all  information  about  the  resolvent  has  already 
been  utilized  in  obtaining  the  characterizer  equation  and  no  informa¬ 
tion  at  all  la  available  for  the  specification  of  the  determinator 
coefficients. 

This,  however,  is  clearly  as  it  should  be,  for  two  unknown 
quantities  were  Introduced  in  Eq.  (8.1)  to  replace  the  one  unknown 
H(u,s;\).  The  determinator  coefficients,  exclusive  of  pQ,  are  therefore, 
in  fact,  entirely  arbitrary.  Any  chosen  set  of  coefficients  pn  may  be 
used  in  Bq.  (8.7)  to  obtain  the  corresponding  set  of  characterizer  coef¬ 
ficients  Cn(u,s)  and,  from  these,  the  resolvent  H(u,s;X.).  Hie  resultant 
expression  will,  of  course,  be  valid  for  only  that  range  of  the  parame¬ 
ter  X  for  which  both  series  of  Eqs.  (8.5)  and  (8.6)  are  convergent. 

In  view  of  the  arbitrariness  of  the  determinator  coefficients, 
it  is  tempting  to  simplify  the  expressions  by  selecting  zero  as  the 
value  of  all  determinator  coefficients,  except  pQ.  This  eliminates  the 
question  of  the  convergence  of  the  series  for  the  determinator  and  sim¬ 
plifies  the  recursive  characterizer  equation  to 

CQ(u,s)  ■  $  K(u,v)Cn_1(v,s)  (8.9) 

which,  with  Co(u,s)  ■  K(u,s),  readily  yields  each  characterizer  coeffi¬ 
cient  in  succession.  With  this  choice  for  the  pQ,  the  final  expression 
for  the  resolvent  is 

H(u,s;\)  -  Jk  K“+1(u,s)Xn  (8.10) 

where  the  "powers"  of  the  kernel  are  defined  recursively  by 

rW)  -  K(u,s)  K3*1^)  -  $  K(u,v)l^(v,s)  (8.11) 
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®iis  is  certainly  a  proper  solution  to  the  resolvent  equation. 
It  is,  in  fact,  the  Neumann  series  solution  which  could  have  been  ob¬ 
tained  directly  from  the  resolvent  equation,  Eq.  (7.3),  by  iteration; 
i.e.,  by  introducing  the  entire  right-hand  side  of  thi3  equation  into 
the  summand  appearing  therein,  and  repeating  the  process  indefinitely. 
This  Neumann  solution  has  the  drawback,  however,  that  the  range  of  X  for 
which  it  converges  is  too  small.  The  series  of  Eq.  (8.10)  actually 
diverges  for  all  values  of  X  with  absolute  value  greater  than  that  of 
the  characteristic  value  of  smallest  magnitude.  As  has  been  remarked, 
the  characteristic  values  of  X  and  the  associated  solutions  of  the  homo¬ 
geneous  equation  are  of  the  greatest  interest,  but  the  Neumann  series 
fails  to  converge  as  soon  as  the  first  characteristic  value  is  attained. 
This  solution  is  thus  entirely  useless  for  the  study  of  the  sourceless 
solutions . 

5Ms  is  the  crux  of  the  problem.  Although  the  determinator 
coefficients  are  arbitrary,  for  a  given  assignment  of  values  to  the  pn 
the  series  for  p(x)  and  for  C(u,s;X)  will  converge  only  for  X  within 
sane  finite  circle  of  convergence  and  the  resulting  expression  for  the 
resolvent  will  be  valid  only  in  some  restricted  range  of  the  parameter 
X.  This  situation  clearly  defeats  the  purpose  of  expressing  the  resol¬ 
vent  in  terms  of  a  characterizer  and  determinator.  What  was  intended 
was  that  the  poles  of  the  resolvent,  which  correspond  to  the  character¬ 
istic  values  of  X  and  the  sourceless  solutions,  be  obtainable  as  the 
roots  of  the  determinator.  The  secular  equation  would  then  be  simply 

p(x)  -  0  (8.12) 

with  an  infinite  number  of  roots  corresponding  precisely  to  the  infinite 
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number  of  characteristic  values  of  X.  associated  with  a  nondegenerate  ker¬ 
nel.  Furthermore,  the  condition  that  Eq.  (8.12)  specify  the  characteris¬ 
tic  values  of  X.  was  to  lead  directly  to  the  actual  nontrivial  sourceless 
solutions  as  well.  Thus,  if  X.Q  solves  Eq.  (8.12),  then  from  Eq.  (8.2) 
C(u,s;\o)  a  \Q  $  K(u,v)C(v,s;\o)  (8.13) 

and  a  solution  to  the  homogeneous  equation,  Eq.  (7.8),  is 

V(u)  -  C(u,s0;\o)co  (8.1k) 

where  cq  is  an  arbitrary,  constant  three-vector  and  sq  is  any  value 
of  the  composite  index  s.  This  sourceless  solution  will  be  nontrivial, 
provided  merely  that  sQ  and  cq  be  not  chosen  so  as  to  annihilate  the 
resulting  V(u). 

To  fulfill  these  desired  conditions,  a  set  of  values  for  the 
determinator  coefficients  must  be  so  specified  that,  firstly,  both  p(\) 
and  C(u,s;X.)  be  entire  functions  of  \,  thereby  insuring  convergence  for 
all  X.  and  the  validity  of  the  expression  for  the  resolvent  over  the  full 
frequency  spectrum;  secondly,  that  the  secular  equation,  Eq.  (8.12), 
yield  all  characteristic  values  and  only  the  characteristic  values  of  X.; 
and  thirdly,  that  the  eigensolutions  C(u,so;X.q)co  be  obtained  as  well. 
How  in  the  degenerate  case,  in  which  p(X.)  is  det(l  -  \  R),  these  condi¬ 
tions  are  satisfied.  To  achieve  these  results  in  the  general  case,  the 
determinator  must  be  specified  in  the  same  way  as  in  the  degenerate  case, 
despite  the  infinite  order  of  the  obverse  matrix.  Wiat  is,  the  relation 
which  is  to  be  specified  so  as  to  determine,  in  conjunction  with  Eq.  (8.2^ 
the  desired  form  of  the  determinator  must  be  precisely  that  which  obtains 
in  the  degenerate  case,  independently  of  the  order  of  the  obverse  matrix. 
This  relation  must  be  investigated. 
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9.  TERMINATION  CONDITION 

By  the  foregoing  reasoning,  the  problem  has  been  reduced  to 
that  of  specifying  a  relation  between  the  character! zer  and  the  deter¬ 
minator  of  the  resolvent  for  a  general  kernel  which  will  yield  the 
optimum  set  of  coefficients  for  the  expansion  of  the  determinator.  This 
optimum  condition  will  be  such  that  the  poles  of  the  resolvent  will  be 
the  roots  of  the  determinator,  both  the  characterizer  and  determinator 
being  entire  functions  of  the  parameter  X.  The  relation  to  be  found 
will  complement  that  of  Eq.  (8.2)  and  thereby  provide  a  sufficient  num¬ 
ber  of  conditions  for  the  unique  determination  of  both  the  characterizer 
and  the  determinator. 

In  the  degenerate  case,  the  determinator  and  the  characterizer 

are 

p(X)  =  det(l  -  X  R)  C(u,s;X)  =  A(u)Q(X)B(s)  (9.1) 

where  Q(X)  is  the  adjoint  of  the  matrix  (1  -  X  R): 

Q(\)  -  (1  -  X  R)_1det(l  -  X  R)  (9.2) 

Clearly,  there  is  an  intimate  relation  between  p(X)  and  C(u,s;X),  as  a 
result  of  which,  effectively,  both  of  these  are  entire  functions  and  the 
characteristic  values  are  given  by  the  roots  of  p(x),  regardless  of  the 
order  of  the  obverse  matrix,  R.  As  expressed  in  Eq.  (9.1),  hcwever> 
both  p(x)  and  C(u,sjX)  depend  on  R  and  can  be  calculated  explicitly  only 
if  its  order,  r,  is  finite.  In  order  to  adapt  the  specification  of  Eq. 
(9.1)  to  the  case  of  a  general  kernel,  it  is  necessary  to  express  the 
relation  between  the  characterizer  and  determinator  in  a  form  which  is 
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independent  of  the  order,  r,  of  the  degeneracy.  Such  a  relation  will  "be 
applicable  to  the  case  of  an  obverse  matrix  of  infinite  order  as  well  as 
to  the  degenerate  case.  The  determinator  and  characterizer  of  a  degen¬ 
erate  kernel  are  entire  functions  because  they  are  polynomials  in  k;  in 
the  nondegenerate  case  their  expansions  in  powers  of  \  will  not  termi¬ 
nate,  but  if  the  condition  on  the  determinator  is  taken  as  the  analog  of 
that  which  holds  in  the  degenerate  case  no  singularities  will  be  intro¬ 
duced  into  either  the  characterizer  or  determinator  to  limit  their  radii 
of  convergence. 

'Hie  relation  sought  is  one  between  the  scalar  p(\)  and  the  me¬ 
tric  function  of  two  composite  indices  C(u,s;k).  The  composite  indices 
can  be  readily  eliminated  by  summing  over  them;  there  will  then  remain 
to  extract  a  scalar  from  the  residual  matrix.  A  significant  set  of 
scalars  associated  with  a  matrix,  closely  related  to  its  eigenvalues  but 
more  easily  calculated,  consists  of  the  traces  of  the  matrix.  An  n  x  n 
matrix  A  is  characterized  by  n  traces.  The  m-th  order  trace,  TrmA, 
is  the  sum  of  the  principal  minors  of  order  m  of  the  matrix.  Its 
significance  lies  in  the  fact  that  TrmA  equals  the  sum  of  the  products 
of  the  eigenvalues  taken  m  at  a  time.  Hie  trace  of  A,  i.e.,  the  sum 
of  its  diagonal  elements,  Tr  A,  corresponds  to  Tr^A;  the  determinant 

det  A  coincides  with  Tr  A.  For  convenience,  the  further  definitions 

n 

will  be  made  that  Tr  A  -  1  and  that  Tr  A  -  0  far  any  m  which  exceeds 

o  m 

the  order  n. 

From  the  relation  between  the  traces  and  eigenvalues  of  a  ma¬ 
trix,  the  following  properties  of  the  traces  may  readily  be  demonstrated. 


Tr  AB  ■ 

Xu 


(9.3) 


Tr  BA 
m 


kh 

Tr  eA  ■  cm  Tr  A  (9.4) 

in  in 

’V1  -  ^n-m  *  /  4«*  *  (9-5) 

f1  +  A>  -  Jo  (nS)  ''k*  (9-6) 

The  conmutatlon  property  of  Bq.  (9*3)  holds  even  If  A  and  B  are  not 
Individually  square  matrices,  although  they  must  be  conformable  in  both 
orders.  This  property  results  in  the  invariance  of  the  traces  under 
similarity  transformations.  In  Bq.  (9.4),  c  is  a  scalar  multiplier. 

Virtually  the  only  significant  scalar  which  can  be  extracted 
from  C(u,s;X)  by  linear  operations  is 

f(\)  -  $  Tr  C(s,s;X)  (9.7) 

It  may  be  expected  that  this  scalar  function  of  X  is  related  to  p(X). 
ibis  relation  constitutes  the  Termination  Theorem,  which  can  now  be 
proved,  granted  the  preliminary  lemma  that  the  operators  $  and  Tr 
commute.  There  follows 

f(x)  -  $  Tr  A(s)Q(X)B(s)  -  $  Tr  Q(x)B(s)a(s) 

-  Tr  $  Q(x)B(s)A(s)  -  TT  Q(X)  R 

where  the  commutation  property  of  Bq.  (9.3)  and  the  linearity  of  $  have 
been  used.  Let 

G(X)  -  -  X  R  (9.8) 

Hence,  by  Bqs.  (9*4),  (9.2),  and  (9.1), 

-  X  f(X)  -  Tr  QQ  -  Tr  (l+G)-1pG  -  pTr  (l+G-1)"1 

Using  Bq.  (9.5)  and  noting  that  det(l-K}“^)  ■  det  [Q"1(l+Q)] 

•  det(l+0)  /  detO  -  p  /  det  0, 
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-  X.  t(\)  .  pTrr_1(l+0-1)/4et(l+0"1)  -  detO  Tr^fl+O-1) 
Proa  Bqs.  (9.6)  and  (9.5), 


r-1 


-  X.  f(X)  -  detG  Z  (r-k)  Tr.  G 


-1 


k-o 

r-1 


-  detG  kZo  (r-k)  Trr_k  G  /  detG 
r 

«  Z,  m  Tr  G 
n*x  m 

which,  from  Biis.  (9*8)  and  (9.k),  Is  equivalent  to 
f(\)  -  -  (-1)“  Trm  R  m  X.®'1 


which  may  he  recognized  to  he 

tM  - 


A*  JL  (*i)ra  r  k' 

dX.  nmo  a 


m 


But,  from  Eq.  (9*6), 


p(X)  -  det(l  -  X.  R)  -  ^  (-1)®  Trm  R  Xm 


Hence  f(\)  -  -  dp/dX.  and  the  relation  sought  is 


$  Tr  C(s,s;\) 


(9.9) 


(9.10) 


which  depends  explicitly  on  neither  R  nor  its  order,  r. 

Conversely,  it  may  he  shown  that  the  power  series  for  hath  the 
determinator  and  characterlzer  will  indeed  terminate  for  a  degenerate 
kernel  if  the  condition  of  Eq.  (9.10)  is  imposed  in  addition  to  Eq. 
(8*2).  For  In  the  degenerate  case,  Bq.  (8.2)  alone  yields  as  the  char¬ 
acterlzer,  for  any  p(\), 

C(u,s;\)  -  A(u)(l  -  XR)"1p(X,)B(s) 
that  $  Tr  C(s,s;\)  -  p(X.)  Tr  (1  -  XB)’1R,  from  which 
-(l/p) (dp/dX)  -  -  (d/dX)log  p  -  Tr(l  -  XR)_1B 
-  -  Tr  (d/dX)log  (l  -  XB) 


(9.H) 


so 


■unm 


or,  in  view  of  the  commutativity  of  Tr  and  d/dX  and  of  the  initial  con¬ 
ditions  on  p(x). 


log  p(X)  -  Tr  log  (1  -  XR) 


(9.12) 


But  since  Tr  log  A  -  log  det  A,  as  follows  from  the  relation  between 
the  trace  and  eigenvalues  of  the  matrix  log  A,  Eg.  (9.12)  implies 

p(x)  -  det  (1  -  XR)  (9.13) 

so  that  the  power  series  for  p(x)  will  terminate  at  Xrj  in  fact. 


p(X)  -  E  (-1)“  Tr  R  Xm 

mao  m 

21  22 

The  Cayley-Hamilton  theorem  *  may  be  written 

E  (-l)r"m  Tr  R  Rm  =  0 
nteo  r‘m 

or,  by  comparing  the  coefficients  in  Eqs.  (9.1U)  and  (9.I5), 

E  p  R®  o  0 


(9.14) 


(9.15) 


By  iterating  in  Eg.  (8.7),  there  is  obtained 
C  (u,s)  «  E  Kmfl(u,s)  p 

X  mao  r' 

But,  from  Eg.  (8.11), 

Kn+1(u,s)  -  A(u)RnB(s) 


(9.16) 


(9.17) 


so  that 


C  (u,s)  -  A(u)  E  p  R®  b(s)  -  0 


(9.18) 


Since  p.^  ■  0  and  Cr(u,s)  ■  0,  Eq.  (8.7)  implies  that  all  higher- 
order  coefficients  of  C(u,s;X)  vanish.  Hence  both  the  determinator  and 
character! ser  expansions  for  a  degenerate  kernel  will  terminate  as  a 


result  of  the  application  of  Eqs.  (8.2)  and  (9.10). 

If  the  termination  condition  13  applied  to  the  nondegenerate 
case,  the  series  vlll  not  actually  terminate,  hut  its  effect,  which  is 
to  cancel  the  poles  of  (1  -  \R)“^  with  the  zeros  of  p(\)  in  Eq.  (9. 11) 
so  as  to  render  C(u,s;X)  an  entire  function,  will  persist  even  as  the 
order  of  the  degeneracy,  r,  is  allowed  to  become  infinite.  Thus,  the 
termination  condition,  Eq.  (9.IO),  realizes  the  proper  pole  cancella¬ 
tion  property  of  the  optimum  determinator. 


10.  COMPLETE  SOLUTION 
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The  charset erizer  equation  and  the  termination  condition, 
together  with  the  initial  conditions,  provide  the  complete  solution  to 
the  problem  posed  herein.  These  equations  are 

C(u,s;X)  ■  K(u,s)p(X)  +  X  $  K(u,v)C(v,s;X)  (10.1) 

-  §£  -  $  Tr  C(s,s;X)  (10.2) 

p(0)  ■  1  C(u,s;0)  *  K(u,s)  (10.3) 

Upon  solving  these  for  the  characterizer  C(u,s;\)  and  determinator  p(X), 
the  resolvent  H(u,s;X)  is  obtained  as  their  ratio,  Eq.  (8.1).  This 
immediately  yields  the  field  transforms  V(s)  and  I(s)  from  Eqs.  (7*2) 
and  (6.4)  and  these  may  he  inverted  to  obtain  the  actual  fields  in  co¬ 
ordinate  space. 

Perhaps  of  greater  importance  than  the  fields  for  the  study  of 
a  given  linear  medium  is  the  fact  that  the  eigenfrequencles  are  obtain¬ 
able  from  the  secular  equation 

p(X)  -  0  (10.4) 

This  equation,  the  central  result  of  the  theory,  has  the  significance  of 
a  dispersion  relation  far  a  propagation  problem,  or  of  a  resonance  condi¬ 
tion  for  a  cavity,  or  of  any  equivalent  relation  among  the  frequency  of  a 
wave,  its  wavelength  and  other  characteristics,  and  the  dimensions  and 
electrical  properties  of  the  medium.  As  such,  the  secular  equation 
affords  a  powerful  characterisation  of  the  electromagnetic  properties  of 
the  medium. 


In  conjunction  with  the  solution  to  the  secular  equation, 
there  is  obtained  in  this  formulation  the  associated  characterlzer, 
which  corresponds  essentially  to  the  residue  of  the  resolvent  at  its 
pole.  This  characterlzer  yields  the  solutions  to  the  homogeneous  equa¬ 
tion,  Eq.  (7.8),  as 

V(s)  -  C(s,so;\o)co  (10.5) 

where  is  a  root  of  the  secular  equation,  sQ  is  any  value  of  the  com¬ 
posite  index,  and  cQ  is  an  arbitrary  constant  three-vector.  The  constant 
vector  co,  together  with  the  selected  value  of  sq,  plays  the  role  of  the 
arbitrary  amplitude  coefficient  for  the  normal  mode. 

The  actual  solution  of  Eqs.  (10.1)  and  (10.2)  may  proceed  by 
recursion.  The  determinator  and  characterlzer  are  expanded  in  power 
series  in  X,  as  in  Eqs,  (8.5)  and  (8.6),  which,  upon  substitution  into 
Eqs.  (10.1)  and  (10.2)  and  comparison  of  coefficients,  yield 


Cn(u,s)  -  K(u,s)pn  +  $  K(u,v)Cn-1(v,s) 

(10.6) 

-npn  -  $  Tr  Cn-1(s,a) 

(10.7) 

P0  -  1  CQ(u,s)  -  K(u,s) 

(10.8) 

These  equations  suffice  to  yield  all  the  unknown  coefficients  of  both 
p(X)  and  C(u,sj\)  in  succession.  Starting  with  Eq.  (10.8),  p^  is  ob¬ 
tainable  from  Eq.  (10. 7),  whereupon  C^(u,s)  can  be  calculated  from  Eq. 
(10.6)  and  this,  in  turn,  will  yield  Pg  through  Eq.  (10. 7);  the  process 
may  be  repeated  indefinitely  to  obtain,  in  principle,  explicit  expan¬ 
sions  of  the  determinator  and  the  characterlzer.  The  rapidity  of  con¬ 
vergence  of  these  expansions  depends,  of  course,  on  the  original  choice 
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of  expansion  functions,  c(r,s),  for  the  transformation.  The  solution  of 
the  problem  would  then  proceed  with  the  solution  and  study  of  the  secu¬ 
lar  equation,  Eq.  (10.2+) . 

There  can  arise  a  situation  in  which  the  procedure  just  de¬ 
scribed  may  be  thwarted  at  its  start.  This  is  the  case  if  K(u,s)  is  in¬ 
finite  or  undefined  for  u  =  s,  whereupon  the  prescription  for  calculat¬ 
ing  p1  from  Eq.  (10. 7)  is  ineffectual.  This  situation  demands  a  modi¬ 
fication  of  the  recursion  relations  which  will  sidestep  the  difficulty. 

That  such  a  modification  is  possible  is  the  result  of  the 
residual  ambiguity  of  the  characterizer  and  determinator.  The  essential 
quantity  is  the  resolvent  and  this  has  been  found  as  the  ratio  H(u,s;\)  * 
C(u,s;\)/p(\).  Both  numerator  and  denominator  of  this  ratio  may  be 
multiplied  by  some  scalar  function  g(\)  to  yield  the  equivalent  ratio 

H(u,s;\)  -  D(u,sj\)/q(\)  (10.9) 

where 

q(\)  -  p(X)g(*.)  D(u,s;\)  -  C(u,s;\)g(\)  (10.10) 

To  preserve  the  initial  value  of  the  determinator,  g(o)  **  1  may  be  speci¬ 
fied.  The  modified  characterizer  equation  is  then  essentially  the  same, 

D(u,sj\)  *  K(u,s)q(\)  +  \  $  K(u,v)D(v,s;\)  (10.11) 

but  the  termination  condition  is  altered  to 

-  ^  $  Tr  D(s,s;\)  -  t(\)q(\)  (10.12) 

where 

t(\)  =  (l/g)(dg/d\)  (10.13) 

Now  if  g(\)  is  chosen  such  that 

t(\)  -  t  •  $  Tr  K(s,s) 


(10.12+) 
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then  the  modified  termination  condition  will  be 

-  -  $  Tr  [D(s,s;X.)  -  K(s,s)q(x)]  (10.15) 

If  power  series  expansions  of  D(u,s;X)  and  q(\)  are  made,  the  modified 
recursive  character! zer  and  termination  equations  become 

Dn(u,s)  ■  K(u,s)qn  +  $  K(u,v)Dn-1(v,s)  (10.16) 

-  n  qn  -  $  Tr  [Djj^Csjs)  -  Kfs,*^^]  (10.17) 

qQ  -  1  Dq(u,s)  •  K(u,s)  (10.18) 

The  result  is  that  ^  •  0  and  manipulation  with  the  possibly  nonex¬ 
istent  K(s,s)  has  been  avoided.  From  Eqs.  (10.13)  and  (10.14)  there 
follows  that 

q(\)  -  etKvM  (10.19) 

so  that  no  new  zeros  have  been  introduced  into  the  determinator,  both 

q(\)  and  D(u,s;\)  are  entire  functions,  and  the  modified  secular  equation 

q(\)  -  0  (10.20) 

yields  the  same  results  as  does  Eq.  (10.4). 

A  more  serious  crisis  can  arise  in  the  previously  described 
method  of  extraction  of  normal  modes,  Eq.  (10. 5).  Ihe  trivial  zero  vec¬ 
tor  may  be  obtained  if  cQ  and  sQ  are  injudiciously  chosen,  but,  normally 
such  values  of  sQ  are  the  exceptional,  easily  avoided  ones.  In  some 
cases,  however,  the  process  may  fail  entirely  far  all  values  of  sQ  and 
cQ  because  C(u,s;\Q)  may  vanish  identically.  This  situation  may  arise 

if  \  is  a  multiple  root  of  the  secular  equation.  The  procedure  for  ex- 
o 

tracting  the  normal  modes  then  becomes  considerably  mare  complicated. 
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Since  no  information  was  lost  in  translating  the  problem  from  Maxwell's 
equations  to  the  present  formulation,  it  may  yet  be  expected  that  some 
further  set  of  operations  may  be  applied  to  the  characterizer  to  extri¬ 
cate  the  complete  set  of  normal  modes  from  the  equations  even  in  case 
is  a  repeated  root  of  p(\). 
purpose  will  now  be  expounded. 


The  appropriate  operations  for  this 
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11.  NOFMAL  MODES 

Die  present  reformulation  of  Maxwell’s  equations  Incorporates 
all  the  results  obtainable  for  any  linear  medium  in  the  corresponding 
resolvent  H(u,s;\).  The  physical  properties  of  the  medium  are  reflected 
in  the  analytic  properties  of  this  function.  In  particular,  the  reso¬ 
nance  condition  or  dispersion  relation  characterizing  the  medium  is  a 
mathematical  statement  which  locates  the  poles  of  the  resolvent,  while 
its  residues  describe  the  normal  modes  of  the  system. 

The  basis  for  obtaining  the  normal  modes  is  the  principle  that 
if  there  be  extracted  from  the  resolvent  any  matrix  T(u,s;X)  with  the 
property  that 

T(u,sj\o)  -  \0  $  K(u,v)T(v,s;\o)  (ll.l) 

where  X  is  a  root  of  the  secular  equation,  then  this  matrix  will  gener- 
o 

ate  solutions  to  the  homogeneous  equation,  Eq.  (7.8),  as  follows. 


V(u) 

-  $  T(u,s;ko)c(s) 

(11.2) 

I(u) 

-  $  Y(u,s)V(e) 

(11.3) 

Here,  c(s)  is  an  arbitrary  three-vector  function  of  one  composite  index. 
Per  this  sourceless  solution  to  be  nontrivial,  it  is  required  merely 
that  the  vector  c(s)  be  not  one  that  is  annihilated  in  sunning  with 
T(u,s;X.Q).  nils  requirement  is  normally  so  weak  that  the  sumaatlon  may 
be  avoided  entirely  by  choosing  a  unit  function  for  the  arbitrary  vector: 
c(s)  ■  col(s,8o),  where  sQ  is  any  value  of  the  composite  index  and  cQ 
any  constant  three-vector,  each  selected  so  that  the  resulting  V(u)  does 
not  vanish.  This  procedure  offers  no  difficulty  unless  T(u,s;\0)  is 


identically  zero. 

As  indicated  previously,  the  characterizer  C(u,sj\)  is  such  a 
matrix  and  it  does  yield  a  normal  mode  for  each  root  of  the  determinator 
in  this  manner.  The  process  fails,  however,  if  C(u,s;\)  vanishes  iden¬ 
tically  at  \  o  \Q.  This  can  occur  if  is  a  multiple  root  of  the 
secular  equation.  Another  matrix  T(u,s;\)  must  then  he  found  which 
satisfies  Eq.  (11.1)  for  this  X.  *  \  and  which  does  not  vanish.  Clearly, 
this  corresponds  to  finding  the  residue  of  the  resolvent  at  a  pole  of 
higher  order.  The  procedure  for  finding  the  normal  modes  for  multiple 
eigenvalues  is,  in  fact,  analogous  to  that  of  extracting  such  a  residue 
and  is  Justified  by  the  following  Successor  Theorem. 


Let  a  differential  operator  D^(x.)  he  first  defined  by 


Dn(\)f(\) 


X.n  dnf 

ST  ix“ 


(11.4) 


It  is  readily  verified  that  this  operator's  effect  on  a  product  is  ex¬ 


pressed  by 

Da(\)  [f(4)g(\)J  -  Jo  VmMf(x)  Dm(\)g00  (11.5) 


Furthermore,  by  virtue  of  the  property  of  the  resolvent  given  in  Eq. 
(7.6),  and  defining  the  "powers"  of  the  resolvent  just  as  those  of  the 
kernel  were  defined  in  Eq.  (8.11),  it  is  found  that 

Dn(\)H(u,s;\)  -  \n  Hn+1(u,sjX)  (U.6) 

Let  new  the  successors  Sn(u,s;\)  and  the  terminators  g^k)  be  defined 

by 

Sn(u,sj\)  -  Dn(\)C(u,s;\)  gn(X)  -  Dn(\)p(\)  (H.7) 

Ihe  relation  between  these  quantities  is  revealed  by  applying  the  Da(\) 
operator  to  the  defining  equation 
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C(u,s;X)  -  H(u,s;X)p(x)  (11.8) 

to  obtain,  by  use  of  Eqs.  (11.5)  and  (U.6), 

Sn(«,s;\)  -  Jo  \n“B  ^•“(u.sj^)  gjx)  (11.9) 

Comparing  this  expression  with  the  same  one  written  for  Sn  ^(u,sjX), 
there  is  immediately  obtained  the  recursion  relation 

Sn(u,sj\)  ■  A.  $  H(u,v;\)Sn  l(v,s;\)  +  H(u,s;X.)gn(\)  (11. 10) 

The  effect  of  the  kernel  K(u,s)  operating  on  a  successor  is  obtainable 
from  this  recursion  relation  by  applying  the  resolvent  equation,  Eq. 

(11.11) 
(11.12) 

Thus,  the  successor  Sn(u,s;\)  will  satisfy  Bq.  (11. l),  provided  that 

R  (u,s;X  V  vanishes, 
n  o 

To  apply  this  Successor  Theorem  to  the  problem  of  calculating 
the  normal  modes  of  a  system,  consider  a  case  of  a  determinator  p(X) 
which  has  a  zero  of  multiplicity  z  at  X  -  XQ;  that  is,  p(X)  has 
(1  -  x/XQ)z  as  a  factor.  It  may  then  be  that  C(u,s;X)  contains 
(1  -  x/XQ)e  as  a  factor,  with  0  <  e  <  z.  Ife/o,  C(u,s;Xo)  vanishes 
and  no  normal  modes  are  obtainable  from  the  characterlzer.  But  the 
appropriate  successor  may  then  be  used  in  lieu  of  the  characterlzer  to 
generate  a  sourceless  solution.  For,  in  this  case,  gn(\Q)  -  0  far 
0  <  n  <  z  while  g„(X  )  i  0  and  S  (u,s;\  )  -  0  for  0  <  n  <  e  while 
Se(u,sjXQ)  ji  0.  This  implies  that  Re(u,s;X0)  ■  0,  so  that  Se(u,sjXQ) 
satisfies  Eq.  (11. l)  and  yields  a  solution  to  the  homogeneous  equation. 


(7.3),  directly.  The  result  is 

Sn(u,sj\)  a  X  $  K(u,v)Sr(v,s;X)  +  Rn(u,s;X) 

Rn(u,s;X)  «  X  $  K(u,v)Sn_1(v,s;X)  +  K(u,s)gQ(x) 
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While  this  disposes  of  the  problem  of  obtaining  a  normal  mode 
at  any  eigenfrequency,  it  leaves  the  question  of  whether  any  other  nor¬ 
mal  modes,  independent  of  the  one  obtained  as  above,  may  exist  for  the 
same  frequency.  The  rather  complex  answer  to  this  question  may  best  be 
derived  from  a  detailed  study  of  the  structure  of  the  resolvent  for  a 
degenerate  kernel,  which  will  now  be  outlined. 

In  the  degenerate  case,  all  sourceless  solutions  are  of  the 

form 

V(s)  =  A(s)  X  (1  -  \  R)  X  >  0  (11.13) 

Hence,  the  question  reduces  to  that  of  determining  the  number  of  indepen¬ 
dent  eigenvectors  X  of  the  obverse  matrix  R  for  \  m  \q.  This  question, 
together  with  the  main  problem  of  how  to  extract  these  normal  modes  from 

the  kernel,  may  be  resolved  by  examining  the  Jordan  canonical  form  of 
21  22 

the  obverse.  *  If  S  is  the  similarity  transformation  matrix  which 
reduces  the  obverse  to  its  canonical  form,  it  becomes  clear  from  an 
examination  of  the  detailed  structure  of  the  matrices  involved  that  the 
normal  modes  sought  are  linear  combinations  of  certain  specific  columns 
of  the  matrix  A(s)S,  the  eigencolumns.  While  this  matrix  is  beyond  cal¬ 
culation,  the  internal  structure  of  the  successors  is  such  that  they  too 
are  formed  from  linear  combinations  of  certain  columns  of  this  matrix. 

The  successors  for  a  degenerate  kernel  are,  from  Eq.  (9.1),  Sn(u,s;\) 

■  A(u)  Dn(\)Q(\)  B(s)  and  upon  examining  the  canonical  form  of  Dn(\)Q(\), 
it  is  found  in  general  that,  first,  for  e  <  n  <  z,  each  sn(u,8;kg)  Is 
formed  from  linear  combinations  of  columns  of  A(u)S  which  1  unhide  some 
eigencolumns;  second,  Se(u,s;ko)  is  composed  of  only  eigencolumns  but 
the  other  successors  include  extraneous  columns;  third,  eigencolumns  not 
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Included  among  the  first  few  nonvanishing  successors  will  he  found  among 
the  remaining  successors;  finally,  all  the  eigencolumns  are  represented 
among  the  complete  set  of  successors  Sn(u,s;\Q)  for  e  <  n  <  z. 

Being  thus  assured  that  all  the  independent  normal  modes  are 
obtainable  from  linear  combinations  of  the  several  successors,  there 
remains  only  to  select  the  appropriate  linear  combinations  to  extract 
only  the  eigencolumns,  filtering  out  the  extraneous  columns,  and  yet 
obtaining  all  the  independent  modes.  Ihe  Successor  Theorem  shows  how 
to  accomplish  this. 

First,  as  already  demonstrated,  at  least  one  normal  mode  may 
be  obtained  from  Se(u,s;\o)  as  V^(u)  *  $  Se(u,s;X.o)c^(s).  If  this  suc¬ 
cessor  incorporates  more  than  one  eigencolumn,  then  these  other  modes 
will  be  obtainable  as  well  by  simply  choosing  other  arbitrary  vectors 

c  (s)  to  form  new,  independent  linear  combinations  of  the  eigencolumns 
m 

in  Se(u,s;\o)  in  the  same  way. 

Next,  once  all  the  independent  modes  within  Se(u,s;\Q)  have 
been  exhausted,  the  eigencolumns  in  Se+^(u,s;\o)  may  be  weeded  out  in  a 
similar  manner,  except  that  the  selection  of  arbitrary  vectors  cm(s) 
must  be  severely  restricted  to  avoid  linear  combinations  in  which  the 
extraneous  columns  appear.  From  Eqs.  (11.11)  and  (11.12),  it  is  clearly 
sufficient  that  these  cm(s)  come  from  just  that  3et  of  vectors  which  was 
excluded  in  the  previous  step;  that  is,  these  new  vectors  must  be  select¬ 
ed  from  among  all  those  which  are  annihilated  by  S  (u,s;\  ),  for  then 
4  R  .  (u,s;\  )c  (s)  vanishes  and  a  new  set  of  normal  modes  arises  from 

firj.  o  m 

$  Se+^(u,s;\Q)cm(s).  More  modes,  in  fact  all  of  them,  will  be  similarly 
obtainable  from  the  other  successors  Sn(u,s;\Q)  by  summing  with  vectors 
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cm(s)  which  ere  annihilated  by  Sn_^(u,s;X.o)  and  not  by  Sn(u,s;h0)  but 
which  are  otherwise  arbitrary..  The  process  must  terminate  no  later  than 
when  Sz_^(u,s;\o)  has  been  forced  to  yield  its  eigencolumns. 

Mich  of  the  vagueness  associated  with  this  process  may  be 
eliminated  through  a  system  of  tabulation  which  traces  the  progress 
towards  the  extraction  of  all  the  normal  modes  from  the  successors,  pro¬ 
vides  upper  bounds  on  the  number  of  modes  which  each  successor  may  be 
expected  to  yield,  and  signals  the  end  of  the  process  as  soon  as  all  the 
modes  have  been  obtained,  thereby  averting  a  fruitless  search  for  non¬ 
existent  modes.  The  tabulation  stems  from  a  classification  of  all  pos¬ 
sible  structures  which  any  matrix  may  have  In  its  Jordan  canonical  form. 
With  the  eigenvalue  of  multiplicity  z,  of  the  obverse  matrix  R  there 
is  associated  a  classifier,  h,  consisting  of  z  positive  Integers,  hm, 
with  the  properties 

A*W  -  * 

\  >  Vi  <u-15> 

Hiere  are  as  many  Independent  eigenvectors  of  R,  or  eigencdLumns  of 
A(u)s,  or  normal  modes,  as  there  are  nonzero  elements  of  h.  To  a  mode 
obtained  frcm  Sn(u,sj\o)  there  corresponds  a  nonzero  element  of  h  of 
value  z-n.  lhese  conditions  Jointly  delimit  the  number  of  steps  re¬ 
quired  in  the  process  of  extraction  of  the  solutions  to  the  homogeneous 
equation  from  the  successors. 

With  the  observation  that  the  terminators  and  successors  can 
be  calculated  directly  from  the  determinator  and  characterlzer  coeffi¬ 
cients,  found  from  Bqs.  (10.6 )-(10.8),  as 


the  foregoing  Is  seen  to  constitute  a  complete  program  for  obtaining  all 
the  normal  modes  of  a  given  linear  medium. 
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12.  GESTURAL  SOLUTIONS 

The  theory  presented  herein  not  only  reformulates  the  elec¬ 
tromagnetics  problem  of  a  linear  medium  in  terms  of  a  new  equation  to 
replace  those  of  Maxwell  but  prescribes  a  step-by-step  procedure  for 
solving  this  resolvent  equation  and  for  extracting  all  the  pertinent 
information  about  the  system  as  veil.  The  information  obtainable  from 
this  formulation  excludes  none  that  could  be  obtained  from  Maxwell's 
equations,  but  it  emphasizes  the  dispersion,  resonance,  or  other  exist¬ 
ence  condition  characteristic  of  the  medium  and  the  complete  set  of  nor¬ 
mal  modes  associated  with  it,  rather  than  the  actual  fields  produced  by 
any  distribution  of  sources  and  the  power  and  energy  carried  by  the 
fields.  The  former  are  actually  the  relations  of  greatest  interest  and 
significance  in  most  situations. 

The  recursive  solution  for  these  quantities  embodied  in  Eqs. 
(K).6)-(10.8)  demonstrates  the  possibility  of  obtaining  a  complete, 
exact  solution  to  the  resolvent  equation.  This  particular  method  of 
solution  may,  however,  be  among  the  least  efficient  of  the  many  which 
are  available.  The  resolvent  equation 

H(u,sj\)  .  m  K(u,s)  +  \  $  K(u,v)H(v,s;\)  (12. l) 

is  the  result  of  central  importance  in  the  reformulation  and  should  be 
considered  on  its  own  merits,  without  any  particular  solution  procedure 
appended  to  it.  A  host  of  analytical  tools  may  be  applied  for  its  solu¬ 
tion,  for  which  purpose  a  compilation  of  the  equation's  main  properties 
and  a  few  approaches  to  its  analysis  should  be  useful. 
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The  Neumann  aeries  solution  to  Eq.  (12.1),  obtained  directly 
by  iteration,  has  been  given  in  Eq.  (8.10).  As  already  noted,  this  ex¬ 
pansion  is  convergent  only  for  small  X  and  cannot  be  used  to  obtain  the 
normal  modes.  It  should  be  useful  for  high-frequency  approximations. 

ip 

Furthermore,  various  analytic  continuation  techniques  may  be  applied 
to  this  series  to  extend  its  validity  to  the  entire  frequency  spectrum. 
The  Fredholm  solution,  Eqs.  (10.1)-(10.3),  has  been  designed  at  the  out¬ 
set  for  validity  at  all  frequencies,  both  the  characterizer  and  deter¬ 
minator  being  entire  functions.  In  this  formulation,  the  existence 
conditions  are  expressed  by  the  vanishing  of  the  determinator. 

For  special  forms  of  the  kernel,  a  complete,  exact,  closed 
solution  is  immediately  obtainable.  One  such  type  of  kernel  is  the 
degenerate  one,  which  was  treated  at  length  in  the  foregoing.  Its  re¬ 
solvent  is  given  exactly  by  Eq.  (7*13)  and  the  existence  conditions  are 
just  those  which  an  eigenvalue  of  the  obverse  matrix  must  satisfy. 
Another  type  of  solvable  kernel,  which  may  be  termed  "ideal,"  is  one  of 
the  form 

K(u,s)  >  N(s)  l(u,s)  (12.2) 

The  unit  function  permits  a  direct  solution  to  the  resolvent  equation. 
The  resolvent  of  the  ideal  kernel  is  found  to  be 

H(u,s;X)  -  [1  -  XN(s)]-1  N(s)  l(u,s)  (12.3) 

and  the  existence  condition  is 

det  [1  -  XN(s)]  —  0  (12.4) 

It  should  be  noted  that  the  composite  index  appears  here  explicitly  in 
the  secular  equation.  linear  combinations  of  degenerate  and  ideal  ker¬ 
nels  also  lead  to  dosed  forms  of  resolvents,  with  the  aid  of  some  other 
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properties  of  resolvents  listed  below. 

The  differential  equation  satisfied  by  the  resolvent  Is 

»  $  H(u,v;\)H(v,s;k)  (12.5) 

With  the  Initial  condition 

H(u,s;0)  *»  K(u,s)  (12.6) 

this  may  he  considered  an  alternate  definition  of  the  resolvent.  More 
basic  still  is  the  functional  equation  for  the  resolvent, 

H(u,s;\)  -  H(u, sj§ )  =  (\-5)  $  H(u,v;\)H(v,s;5)  (12.7) 

which  identifies  the  class  of  functions  to  which  all  resolvents  belong. 

General  properties  of  resolvents  include  that  of  conanutation: 

$  H(u,v;X)H(v,s;§)  »  $  H(u,v;§  )H(v,s;\)  (12.8) 

In  particular,  the  resolvent  comnutes  with  the  kernel.  The  resolvent  of 
a  kernel  which  is  proportional  to  one  whose  resolvent  is  known  is  found 
as 

K(u,s)  ■  a  Kq(u,s)  H(u,s;X)  •  a  H0(u,s;a\)  (12.9) 

If  two  kernels  are  orthogonal;  that  is,  if 

*  KjtavjK^s)  »  $  Kg^vjKjKs)  .  0  (12.10) 

then  the  resolvent  of  the  sum  is  the  sum  of  the  resolvents: 

K(u,s)  *  Kj^u^J+Kg^s)  H(u,s;\)  ■  H1(u,s;\)+H^(u,s;\)  (12.11) 

The  resolvent  of  a  kernel  which  is  the  sum  of  a  degenerate  one  and  one 
whose  resolvent  is  known  is  expressible  in  closed  farm,  as  follows.  If 
K(u,s)  -  A(u)B(s)  +  K0(u,s)  (12.12) 

then  let 


Aq(u;\)  -  A(u)  +  \  $  Ho(u,s;\)A(s) 


(12.13) 
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Bq(sj\)  -  B(s)  +  X  $  B(u)H0(u,s;\)  (12.1h) 

*0M  -  $  B(8)Ao(s;\)  -  $  B0(sj\)A(s)  (12.15) 

and  the  resolvent  for  this  partially  degenerate  kernel  is 

H(u,s;X)  -  H0(u,s;X)  +  Ao(u;X)[l-XRo(x)]"1Bo(s;X)  (12.16) 

Finally,  it  is  noteworthy  that  the  recursion  relations  of  Bis. 
(10.6)-(10.8),  which  produce  the  power  series  expansions  for  the  charac- 
terizer  and  determinator  of  the  Fredholm  solution,  are  explicitly  solv¬ 
able,  yielding  expressions  for  each  of  the  determinator  and  characterlzer 
coefficients  which  may  be  evaluated  directly  from  the  kernel.  The  kernel 
must  first  be  processed  to  yield  the  terminants, 

tn  -  (-l/n)  $  Hr  K“( s,s)  (12.17) 

In  terms  of  these,  the  explicit  expressions  for  the  determinator  and 
characterlzer  coefficients  are 

p  .  r  ft  Vi  (12.18) 

B  »- 1  V 

=„(“.»)  -  Jo  k”1  (12-19) 

In  Eq.  (12.18),  the  summation  is  over  all  terms  of  the  indicated  type 
which  may  be  farmed  from  n  positive  integers  e„  such  that  0  <  e_  <  n 
and 

Jloem  "  n  (12.20) 

Eqs.  (12.17)  and  (12.18)  epitomize  the  entire  processing  required  to 
extract  the  existence  condition  from  the  kernel.  The  recursion  rela¬ 
tions  may,  however,  be  more  convenient  than  these  cumbersome  explicit 
expressions. 


3&e  foregoing  brief  summary  of  the  analytic  properties  of 
resolvents  is  intended  as  a  guide  to  the  most  efficient  way  of  attack¬ 
ing  the  resolvent  equation  for  any  particular  problem.  It  may  be  seen 
that  the  ease  with  which  the  equation  may  be  manipulated  and  solved 
depends  mainly  on  the  fcrm  of  the  kernel.  A  crucial  point  to  be  noted 
is  that  the  form  of  the  kernel  is  in  great  measure  dictated  by  the  orig¬ 
inal  choice  of  the  basis  of  representation  for  the  solution.  Although 
the  choice  of  the  complete  sets  is  in  principle  arbitrary,  a  Judicious 
initial  choice  will  clearly  simplify  all  further  calculations. 
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13.  GENERAL  TECHNIQUES 

The  theory  developed  herein  nay  be  extended  through  general 
and  special  techniques  applicable  in  various  situations.  Specialization 
may  be  desirable  in  certain  cases  with  simplifying  features  for  which 
the  use  of  the  formidable  machinery  set  up  here  would  be  an  extravagance. 
Generalizations,  on  the  other  hand,  can  lead  to  the  relaxation  of  some 
of  the  restrictions  under  which  the  theory  is  valid. 

One  restriction  which  may  easily  be  dropped  is  that  which  re¬ 
quires  the  constitutive  tensors  to  be  independent  of  frequency.  Ihis 
condition  was  imposed  for  convenience  of  interpretation  of  the  secular 
equation  as  a  dispersion  relation  or  resonance  condition.  Many  of  the 
equations  obtained,  particularly  the  termination  condition,  are  not 
valid  if  the  kernel  is  a  function  of  the  parameter  which  is  the  case 
when  the  constitutive  tensors  depend  on  the  excitation  frequency. 

To  adapt  the  theory  for  the  case  of  frequency-dependent  con¬ 
stitutive  parameters  requires  only  the  redefinition  of  the  parameter  \ 
and  of  the  kernel  K(u,s).  Abandoning  the  original  relation  between  \ 
and  the  frequency,  expressed  in  Eq.  (6.5)>  the  parameter  X.  is  now  to  be 
considered  an  arbitrary  auxiliary  variable  unrelated  to  any  physical 
quantity.  The  kernel  K(u,s)  of  Eq.  (6.6)  is  new  more  conveniently  re¬ 
defined  as  simply 

K(u,s)  -  $  Z(u,v)Y(v,s)  (15*1) 

The  entire  subsequent  theory  remains  valid,  with  only  the  following 
slight  modification.  With  this  new  definition  of  the  kernel,  Eq.  (6.8) 


66 


represents  a  problem  somewhat  more  general  than  required,  but  one  which 
reduces  to  that  of  Maxwell's  equations  fcr  X.  ■  1.  This  last  specifica¬ 
tion  is  therefore  to  be  Introduced  into  all  the  equations  so  as  to  make 
them  valid  for  the  electromagnetics  problem  in  question..  In  particular, 
the  new  secular  equation  is 

p(l)  -  0  (13.2) 

which,  despite  the  appearance  of  overdetermination,  simply  fixes  the 
wave  number  k  in  relation  to  the  dimensions  and  other  parameters  of  the 
system.  Eq.  (13.2)  can  hence  be  interpreted  in  the  same  way  as  was  Eq. 
(10. k).  Complications  may  arise,  however,  since  the  expansions  of  the 
determinator  and  characterizer  can  no  longer  be  claimed  to  converge  over 
the  entire  frequency  spectrum,  for  the  frequency  dependence  of  the  con¬ 
stitutive  parameters  may  introduce  singularities  into  these  functions. 
Analytic  continuation  techniques  may  then  be  required  to  suit  the  partic¬ 
ular  frequency  dependence  of  the  medium. 

A  more  sweeping  alteration  of  the  theory  is  called  for  in  case 
the  kernel  is  singular.  A  weak  type  of  singularity  has  already  been 
disposed  of  through  a  modification  of  the  termination  condition  as  in 
Eq.  (10.15).  For  stronger  singularities,  however,  in  which  higher  powers 
of  the  kernel  are  also  singular,  this  modification  is  ineffectual.  This 
is  the  case,  for  example,  with  the  ideal  kernel  of  Bq.  (12.2).  What  is 
required  in  such  cases  is  a  redefinition  of  the  determinator,  as  follows. 

H(u,s;\)  -  C(u,s;\)/p(s;\)  (13.3) 

The  ideal  determinator,  p(s;\),  is  new  a  function  of  the  composite  index 
s  as  well  as  of  the  parameter  \.  The  termination  condition  must  corres¬ 
pondingly  be  modified  to 
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“  -aj2kJ?^  -  $u  C(u,sj\)  (13.10 

where,  as  Indicated,  the  generalized  summation  is  to  be  performed  only 
over  the  composite  index  u.  A  similar  modification  in  summation  is  to 
be  introduced  into  the  appropriate  equations  of  the  rest  of  the  theory, 
which  thereupon  remains  valid  even  in  such  singular  cases. 

The  new  secular  equation  for  such  cases  is 

P(sj\)  -  0  (13.5) 

Its  interpretation  is  no  more  difficult  than  that  of  the  previous  form. 
Here,  the  existence  condition  relates  the  wave  number  k  implicit  in  X.  to 
not  only  the  dimensions  of  the  system  but  to  the  transform  index  s  as 
well.  This  index  will  then  readily  be  assigned  a  physical  significance, 
such  as  that  of  a  propagation  vector  in  a  dispersion  relation. 

Specialization  of  the  theory  is  called  for  to  take  advantage 
of  such  simplifying  features  of  a  problem  as  its  symmetry  properties. 

For  example,  a  typical  propagation  problem  may  involve  axial  symmetry, 
which  invites  the  introduction  of  a  propagation  factor  of  the  form 
g-Jfcz  in  cylindrical  coordinates.  This  simplifies  the  subsequent  cal¬ 
culations  since  this  will  be  the  dominant,  if  not  the  only,  axial  vari¬ 
ation  of  the  fields.  The  resulting  simplification  can  be  introduced  at 
the  outset  by  modifying  the  curl  operator  of  Eq.  (3.2)  to  the  extent  of 
replacing  &/&z  by  3/3z  -  jb  and  then  eliding  the  propagation  factor. 
Similar  modifications  may  be  made  in  case  of  circular  or  other  symmetry. 

More  generally,  the  presence  of  symmetry  in  a  problem  can  lead, 
with  the  appropriate  choice  of  the  basis  of  representation,  to  a  reduction 
in  the  dimensionality  of  the  resolvent  equation.  If  the  composite  index 
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can  be  partitioned  as  s  -  (a^,Sg),  so  that  the  kernel  is  ideal  with 
respect  to  one  part  of  the  index,  then  the  multiple  summations  can  be 
reduced  in  number.  Thus,  if  K(u,s)  »  K(s1;u2,8g)l(u1,s1),  then 
H(u,  s;\)  ■  H(s1;Ug,SgjX.)l(u^,s1),  with  the  dimensionally  reduced 
resolvent  equation 

H(si;ug,s2;\)  a  K(s1JUg,Sg)  +  X.  K^jUg^gjH^jVg^gjX.)  (13.6) 

in  which  the  status  of  is  that  of  a  parameter,  rather  than  of  a  dummy 
index. 

Another  feature  of  some  media  which  allows  the  calculations  to 
be  reduced  is  the  presence  of  conducting  regions  or  boundaries.  These 
can  be  treated  formally  by  introducing  infinite  values  of  the  conductiv¬ 
ity  into  the  permittivity  dyadic  for  those  values  of  the  position  vector 
which  correspond  to  the  perfectly  conducting  regions  of  the  space.  In 
the  case  of  a  closed  system  bounded  by  conducting  walls,  the  result  is 
that  the  integrations  involved  in  the  transformations  become  limited  to 
the  interior  region  of  the  system. 

If  approximate  solutions  suffice  for  the  application  at  hand, 
many  simplifying  techniques  are  available.  A  general  kernel  could  be 
approximated  by  one  which  is  degenerate,  whereupon  an  exact,  closed- 
form  solution  for  the  resolvent  for  the  approximate  kernel  is  obtainable. 
Another  situation  in  which  the  full  machinery  of  the  theory  is  needlessly 
cumbersome  if  approximate  solutions  are  adequate  is  one  which  represents 
a  slight  perturbation  of  a  medium  for  which  the  kernel  is  degenerate. 

The  kernel  will  then  be  expressible  as  in  Eq.  (12.12),  in  which  KQ(u,s) 
is,  in  seme  appropriate  sense,  small.  The  Neumann  series,  Eq.  (8.10), 
is  then  a  rapidly  convergent  expression  for  the  resolvent  H0(u,sj\)  for 
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the  perturbing  kernel  Kq(u,s).  The  solution  given  In  Eq.  (12.16)  can 
then  be  used,  with  Hq(u,s;\)  known  at  least  approximately.  Still 
another  type  of  ready-made  approximation  is  available  in  case  the  kernel 
la  sharply  peaked  near  u  ■  s,  though  perhaps  not  actually  of  the  ideal 
type.  The  approximate  resolvent  is  then  easily  seen  to  be 

H(u,sj\)  -  [1  -  k  KtavjrtKs)  (13.7) 

This  last  approach  can  be  the  basis  for  a  method  of  "moments"  far  approx¬ 
imating  the  resolvent,  in  which  the  resolvent  is  expanded  in  a  power 
series  about  the  value  of  the  composite  index  at  which  the  kernel  is 
peaked  and  a  set  of  simultaneous  metric  equations  is  solved  for  the 
resolvent . 

An  extensive  generalization  of  the  theory  may  be  undertaken 
with  the  goal  of  lifting  the  restriction  to  time-invariant  media.  While 
slow  time  variations  of  the  constitutive  parameters  present  no  diffi¬ 
culties  in  yielding  to  analysis  as  modulated  waves,  rates  of  variation 
comparable  to  those  of  the  fields  would  require  the  abandonment  of  the 
original  steady  state,  harmonic  analysis.  The  theory  could  be  reformu¬ 
lated  in  terms  of  a  maxi  operator  which  incorporates  the  partial  deriv¬ 
atives  with  respect  to  time.  The  only  remaining  restriction  upon  the 
applicability  of  the  theory  would  then  be  simply  to  media  which  are 
linear. 
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14.  ILLUSTRATION 

As  an  illustration  of  the  applicability  and  feasibility  of  the 
method  of  analysis  presented  herein,  there  will  now  be  obtained  there¬ 
with  the  solution  to  a  particular  class  of  problems  in  electromagnetic 
theory.  Although  this  example  is  computationally  quite  trivial  and  is 
certainly  not  intractable  to  standard  methods  of  analysis,  it  does  aid 
in  concretizing  the  rather  abstract  theory  and  possesses  sufficient 
generality  to  be  of  interest.  A  specific  example  of  this  class  of  prob¬ 
lems  will  also  be  examined  in  detail  in  order  todLsplay  the  nature  of 
the  calculations  involved  in  the  matrix  formulation  which  Tinder  lies  the 
theory. 

The  system  to  be  analyzed  is  that  of  an  infinite,  homogeneous, 
linear  medium.  Such  a  medium  is  characterized  by  permittivity  and  per¬ 
meability  tensors  which  are  independent  of  position,  though  otherwise 
arbitrary: 

e(r)  =  e  u(r)  *»  u  (l4.l) 

The  complete  solution,  as  embodied  in  the  resolvent  for  this  system, 
will  be  obtained  and,  in  particular,  the  dispersion  relation  for  the 
medium.  The  theory  may  be  applied  as  follows. 

As  there  is  no  characteristic  physical  dimension  of  the  system 
and  since,  moreover,  the  constitutive  tensors  may  be  functions  of  fre¬ 
quency,  the  kernel  as  defined  in  Eq.  (13.1)  will  be  used  and  X.  will  be 
set  equal  to  unity  in  the  final  result.  For  lack  of  any  reason  to  com¬ 
plicate  the  calculations,  the  transformer  c(r,s)  which  farms  the  basis 
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of  representation  will  be  assumed  to  be  a  scalar.  The  lsmlttanoe  ker¬ 
nels  are  then  given  by 

Z(u,s)  ■  $  d(u,r)(ke)**^curl  c(r,s)  (lit-. 2) 

Y(u,s)  ■  $  d(u>r)(k(i)"3‘crarl  e(r,s)  (14.3) 

as  Is  clear  from  an  Inspection  of  Eqs.  (6.2),  (It. Ik),  (3>7)>  (3*4),  and 
(3.5).  As  a  result  of  the  constancy  of  the  constitutive  tensors  and  the 
scalar  character  of  d(u,r),  these  may  be  rewritten  as 

Z(u,s)  ■  (ke)”1  $  d(u,r)  curl  c(r,s)  (14.4) 

Y(u,s)  »  (ku)"1  $  d(u,r)  curl  c(r,s)  (14.5) 

In  view  of  the  orthonormality  relation  of  Eq.  (It. 7),  It  Is  clearly  ad¬ 
vantageous  to  select  c(r,s)  so  that  the  curl  operation  leaves  It  Intact. 
If  rectangular  coordinates  are  selected,  this  is  easily  accomplished, 
far  then  the  operator  of  Eq.  (3.1)  will  clearly  leave  an  exponential 
function  unchanged  in  form.  Hence 

c(r,s)  ■  e“^r*s  (14.6) 

Is  selected,  the  composite  Index  s  being  considered  a  vector.  This 
amounts  to  nothing  more  than  choosing  a  Fourier  transform  to  effect  the 
solution.  The  result  of  operating  on  c(r,s)  with  the  curl  operator  is 


expressed  by 

curl  c(r,s)  -  c(r,s)(-j/) 

(14.7) 

as  follows  from  Eqs.  (3*1)  and  (3*7)*  Hence, 

$  d(u,r)  curl  c(r,s)  -  l(u,s)(-J^) 

(14.8) 

and  the  lsnlttance  kernels  are 

Z(u,s)  -  (-JA)*"1^ 

(14.9) 

Y(u,s)  - 

(14.10) 
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The  kernel  Is  therefore  simply 

K(u,s)  .  (-l/k2^'1^"1/*  l(u,s)  (14. ll) 

JMs  kernel  is  of  the  ideal  type,  as  in  Eq.  (12.2),  with  N(s) 

2  •■il  «»X 

■  (-l/k  )g"  /n’  The  resolvent  for  such  kernels  is  given  in  Eq.  (12.3). 
In  this  case,  letting 

D(s)  -  (14.12) 

the  resolvent  is 

H(u,s;X)  -  [1  +  (\/k2)D(s)]-1(-l/k2)D(s)  l(u,s) 

or,  since  X  ■  1, 

H(u,sjl)  =  -  [k2  +  D(s)]"1D(s)  l(u,s)  (14.13) 

and  the  problem  is  solved. 

The  dispersion  relation  for  this  medium  is  given  by  Eq.  (12.4), 
which  may  be  written  in  this  case 

det  [k2  +  D(s)]  -  0  (14.14) 

Alternatively,  the  ideal  determinator  p(s;X)  may  be  obtained  through  the 
recursion  relations  as  an  expansion  in  powers  of  X,  which  is  unity. 

Vhat  is  obtained  then  is  the  same  as  the  expansion  of  Eq.  (14.14)  by  use 
of  Eqs.  (9.6)  and  (9.4).  The  result  is 

k6  +  TrjDfs)  k4  +  Tr2D(s)  k2  +  Tr^(s)  -  0 

or,  since  Tr^D(s)  =  det  D(s)  =  0  as  a  result  of  the  fact  that  det  /(  *  0, 
k4  +  Tr  D(s)  k2  +  Tr2D(s)  -  0  (14.15) 

This  is  the  dispersion  relation  for  any  infinite  homogeneous  medium. 

From  Eq.  (14.6),  s  is  clearly  the  propagation  vector  and  the  transforma¬ 
tion  is  merely  an  expansion  in  plane  waves.  The  normal  modes  are  given 
by  the  condition  that  V(s)  be  an  eigenvector  of  the  matrix  D(s). 


Dlls  disposes  of  all  Infinite,  homogeneous,  linear  media,  he 
they  anisotropic,  dispersive,  lossy,  with  or  without  electric  and  mag¬ 
netic  current  source  distributions  of  any  kind. 

A  specific  example  of  a  medium  of  this  type  which  is  of  some 
interest  will  illustrate  the  detailed  matrix  calculations  required  by 
this  formulation.  Consider  a  medium  which  is  gyrotropic,  both  in  capac¬ 
itivity  and  in  permeability.  The  former  type  is  realizable  in  a  plasma, 
the  latter  in  a  ferrite.  If  the  preferred  direction  for  both  is  the 
z-axis,  the  constitutive  tensors  take  the  form 


u 


Ml  JMg 

-JM2  Mi 

0  0 


(14.16) 


Hie  explicit  dispersion  relation  for  this  medium  is  desired.  Prom  Eq. 
(14.15),  this  requires  the  calculation  of  the  traces  of  D(s)  - 
A  direct  calculation  involves  considerable  algebraic  manipulation,  much 
of  which  is  superfluous.  By  taking  advantage  of  the  structure  of  the 
constitutive  tensors  and  of  the  matrix  the  process  is  reducible  to 
elementary  matrix  manipulations,  with  a  spectacular  saving  in  labor,  as 
follows. 


The  structure  of  the  constitutive  dyadics  invites  a  partition 
of  the  matrices  into  transverse  and  longitudinal  parts: 


where  Z  and  7  are  2x2  matrices.  Defining  the  numerical  matrix  S  as 


and.  noting  the  properties 


-1 


-S 


(14.19) 


where  +  here  denotes  the  transpose  of  the  matrix,  the  Inner  structure  of 
the  partitions  can  he  conveniently  exhibited: 

Z  1  -  e1-tJCgS  z"1  -  Y’1  -  p^-JpgS  y"1  -  p^  (14.20) 

It  follows  from  Eq..  (14.19)  that  the  inverses  partition  with  the  same 
structure: 

Z  -  z1  +  ZgS  Y  «  yx  +  ygS  (14.21) 

from  which  it  is  clear  that  Z  and  Y  commute  with  each  other  and  with  S. 
This  simple  observation  in  Itself  leads  to  a  considerable  saving  in 
labor. 


Hie  propagation  vector  s  can  also  be  partitioned  into  trans¬ 
verse  and  longitudinal  parts  so  that  the  structure  of  the  ^  matrix  may 
also  be  revealed! 


(14.22) 


where  b  is  a  2  x  1  column  and  c  is  a  single  element.  Hence,  D(s)  is  to 
be  obtained  as 


D(s) 


(14.25) 


which  will  have  the  structure 
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D(s)  . 


(14.24) 


Hie  traces  of  this  matrix  are  expressible  in  terms  of  the  par¬ 
titions  as 

Tr  D(s)  -  d  +  Tr  Dq  (14.25) 

Trg  D(s)  ■  det  Dq  +  d  Tr  Dq  -  dgd^  (14.26) 


Even  the  calculation  of  det  Dq  in  Eq.  (14.26 )  can  be  avoided  by  recall¬ 
ing  that  det  D(s)  is  known  to  be  zero: 

det  D(s)  -  0  -  d  det  Dq  +  dgSD^Sc^  (14.27) 

Hence,  Trg  D(s)  can  be  expressed  as 

•Erg  D(s)  -  (-1/d)  dgSD^Sd1+dTrDo-  d^  (14.28) 

which  involves  nothing  more  difficult  than  matrix  multiplication. 

Now  the  partitions  of  D(s)  are  immediately  obtainable  from 
Eq.  (14.23): 


D  <■ 
o 

c^ZSYS  +  yZSbb+S  -  -c2ZY  +  y!3Zbb+S 

(14.29) 

d  - 

zb+SYSb 

-  -zb+Yb 

(14.J0) 

*1  m 

-cZSYSb 

»  cZYb 

(14.31) 

*2  - 

-zcb+SYS 

-  zcb+Y 

(14.52) 

where  the  commutation  relations  noted  above  and  the  properties  of  Bq. 
(14.19)  have  been  used  to  effect  the  simplifications.  The  slight  labor 
involved  in  evaluating  Eqs.  ( 14.25)  and  (14.28)  can  be  reduced  still 
further  by  noting  that  b+Sb  «  0  since  S  is  antisymmetric,  that  Tr  S  •»  0 
and  Tr  1  >  2,  that  the  trace  is  a  linear  operator,  and  that  matrices 
may  be  commuted  in  taking  traces,  as  in  Bq.  (9*5)*  Making  use  of  Bq. 
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(14.21),  there  is  obtained,  with  little  effort. 


Tr  D(s)  -  2c2(z2y2-z1y1)  -  b^zy^z^ 

0A.33) 

Tr2  D(s)  -  c\y24y2)  (z2+z|)  +  b^zyz^  + 

c^2  [yy1(z2+z2)  +  zz1(y24y2)] 

(14.34) 

2  + 

In  these  expressions,  b  =  b  b  is  the  square  of  the  magnitude  of  the 
transverse  part  of  the  propagation  vector.  The  traces  can  easily  be 
rewritten  in  terms  of  the  elements  of  the  constitutive  tensors  by  noting 
that  z£+Zg  *  det  Z  «  l/det  Z~  ■  l/(«*  -  €g)>  with  a  similar  expression 
involving  Y. 

Introducing  these  traces  into  Eq.  (14.15)  gives  the  dispersion 
relation  for  the  doubly  gyrotropic  medium.  Its  physical  significance 
becomes  more  apparent  when  it  is  written  in  terms  of  the  magnitude,  p, 
of  the  propagation  vector  s  and  of  the  angle,  0,  between  the  direction 

of  propagation  and  the  preferred  direction  in  space. 

.2  .2.2  2  .2  2.  .2  J2  .  2.  ,, , 

p  c  b  +  c  c=3  cos  ©  b  ■>  p  sin  0  (14.35) 

In  terms  of  a  refractive  index  n  -  p/k,  the  dispersion  relation  can  be 
written  after  a  little  rearrangement  as 

An1*  -  Bn2  +  C  -  0  (14.36) 

with  the  coefficients  given  by 

2  2  2  2 

A  a  (e^sin  9  +  e^cos  ©)  (^sin  ©  +  ^cos  ©)  (14.37) 

2 

B  ■  2*^(6^  +  e2Jig)cos  ©  + 

[^^(e^-eg)  +  sin2©  (14.38) 


(lfc.39) 
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Eq.  (14.36)  Is  a  quadratic  in  n  ,  which  may  readily  he  solved  to  yield 
the  refractive  index  in  any  direction. 

Ihe  homogeneity  of  the  space  makes  this  illustration  a  trivial 
application  of  the  theory,  the  full  power  and  generality  of  which  is  only 
barely  in  evidence  in  this  example.  Here  the  integrations  were  so  ele¬ 
mentary  as  to  reduce  the  calculations  to  little  more  than  simple  algebra. 
In  general  inhomogeneous  media  the  quadratures  may  be  inordinately  numer¬ 
ous  and  cumbersome  but  the  mechanics  of  setting  them  up  and  combining 
them  to  extract  the  desired  results  are  much  the  same  as  in  this  simple 


illustration 
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15.  CQRCLU8I0HS 

Hie  formalism  developed  In  this  work  attempts  to  provide  a  new 
vantage  point  from  which  to  view  and  attack  a  large  class  of  problems  in 
electromagnetic  theory.  Concentrating  on  the  invariable,  characteristic 
features  of  general  linear  media  and  subordinating  their  less  signifi¬ 
cant  details,  it  unifies  the  approach  to  a  wide  variety  of  such  media, 
to  which  the  various  standard,  classical  modes  of  attack  have  assigned  a 
variety  of  epithets  —  inhomogeneous,  anisotropic,  etc.  Ihe  method  of 
analysis  reformulates  Maxwell's  equations,  provides  an  alternate  start¬ 
ing  point  for  the  extraction  of  all  significant  information  from  a  given 
system.  Recognizing  that  a  full  description  of  the  medium  must  implic¬ 
itly  contain  all  the  relevant  information  about  the  system,  the  formula¬ 
tion  seeks  a  direct  route  from  an  initial  mathematical  description  of 
the  electrical  constitution  of  the  space  to  all  the  information  about 
the  consequent  phenomena  which  may  be  of  interest  and  significance. 

The  present  reformulation  combines  many  disparate  analytical 
techniques  of  great  power  —  those  of  partial  differential  equations,  of 
matrix  algebra  and  calculus,  of  abstract  linear  operators,  of  general¬ 
ized  transforms,  and  of  Integral  equations  —  which  may  be  brought  to 
bear  upon  the  general  problem.  Acknowledging  the  self-defeating  fea¬ 
tures  of  seeking  closed-form  solutions  for  media  with  any  but  the  sim¬ 
plest  structures,  the  formalism  allows  considerable  freedom  of  choice  of 
the  form  in  which  the  results  are  represented.  Care  has  been  taken  to 
preserve  intact  the  informational  content  of  Maxwell's  equations  in- 
effecting  the  translation  to  the  new  language. 
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Besides  providing  a  single  new  equation  of  radically  different 
farm  to  be  analyzed  for  the  solution  of  a  problem,  this  work  presents  a 
host  of  mathematical  properties  of  this  equation  which  may  be  used  to 
effect  its  solution  or  to  derive  qualitative  features  of  the  results. 

The  method  is  crowned  by  the  prescription  of  an  explicit  sequence  of 
mathematical  operations  which,  in  principle,  lead  ultimately  to  the  com¬ 
plete  solution  for  any  given  linear  medium.  This  formal  solution  in¬ 
cludes  the  fundamental  existence  condition  characteristic  of  the  medium, 
the  complete  set  of  normal  modes  for  the  system,  the  entire  electromag¬ 
netic  field  pattern  in  response  to  any  given  distribution  of  electric 
and  magnetic  current  sources,  and  the  power  flow  and  energy  distribution 
accompanying  the  fields  —  in  short,  all  the  information  which  is  latent 
in  the  description  of  the  medium  and  excitations.  In  effect,  there  is 
set  up  a  fictive  machine  which  incorporates  the  physical  mechanism  of 
Interaction  of  electromagnetic  fields  and  sources,  has  as  its  input  a 
full  pointwise  description  of  the  electrical  content  of  the  space,  and 
produces  as  its  output  all  the  desired  results  describing  the  phenomena 
associated  with  the  medium.  Effectively,  the  formal  solution  unifies 
and  systematizes  the  analysis  of  a  general  class  of  problems  in  electro¬ 
magnetic  theory  and  reduces  such  problems  to  quadratures. 

The  objection  may  be  raised,  and  a  formidable  one  it  may  well 
be,  that  the  quadratures  involved  in  the  process  are  of  such  complexity 
and  multiplicity asto be  prohibitive.  In  the  same  vein,  the  criticism  may 
he  made  that  the  convergence  of  the  expansions  for  the  deterndnator  and 
characterlzer  may  he  so  slow  as  to  require  an  excessive  number  of  terms 
in  the  series,  with  all  the  calculations  which  that  entails.  These 
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objections  may  be  palliated  by  noting  that  the  series  expansion  method 
is  only  one  of  many  modes  of  attacking  the  resolvent  equation  and  is 
generally  not  the  most  efficient  one.  Furthermore,  the  complexity  of 
the  computations  depends  on  the  choice  of  the  basis  of  representation 
for  the  results.  All  possible  a  priori  knowledge  of  the  form  of  the 
solution  should  hence  be  brought  to  bear  upon  the  choice  of  basis  func¬ 
tions.  If  some  members  of  the  complete  set  of  expansion  functions  resem¬ 
ble  the  functions  to  be  expressed,  the  number  of  significant  terms  re¬ 
quired  in  the  expansions  will  be  small,  though  possibly  at  the  cost  of 
complicating  the  integrations  and  summations  involved.  In  practice, 
some  compromise  wild  be  arrived  at  between  rapidity  of  convergence  and 
complexity  of  calculations.  The  dominant  role  which  can  be  played  by 
automatic  computers  in  this  connection  should  be  noted. 

Elaborations  and  refinements  of  the  theory  may  take  several 
directions.  The  formalism  could  be  expanded  to  handle  rapidly  time- 
varying  and  nonstationary  media  in  a  unified  way.  In  fact,  a  relativ- 
istlcally  invariant  formulation  could  be  developed  by  operating  on  four- 
tensors  instead  of  three-vectors.  Although  the  equations  governing  such 
cases  are  of  a  form  quite  different  from  that  of  the  Maxwell  equations 
treated  here,  the  basic  ideas  used  in  this  work  are  evidently  applicable 
to  any  set  of  linear  partial  differential  equations,  whatever  their 
domain  of  definition.  Toward  the  improvement  of  the  efficiency  of  the 
method,  studies  could  be  undertaken  of  means  of  estimating  the  remain¬ 
ders  of  the  series  appearing  in  the  resolvent  when  these  are  truncated, 
of  a  decision  procedure  far  optimizing  the  choice  of  the  basis  of 
representation,  and  of  general  methods  of  improving  the  convergence  of 
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the  series  representation  of  the  results. 

Valuable  applications  of  the  theory  may  be  made  in  the  analy¬ 
sis  of  various  inhomogeneous  and  anisotropic  media  such  as  that  of  a 
plasma.  Ultimately,  this  work  may  lead  to  successful  attacks  upon  such 
fundamental  problems  as  the  formulation  of  the  conditions  for  the  appear¬ 
ance  of  certain  vave  phenomena  and  the  development  of  synthesis  proce¬ 
dures  for  the  implementation  of  desired  dispersion  relations. 
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